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A  branch-and-bound  solution  algorithm  and  a  computer  program 
implementing  this  algorithm  are  developed  to  solve  a  multiactivity 
multifacility  capacity  constrained  0-1  assignment  problem.  The  math¬ 
ematical  formulation  for  such  a  problem,  called  problem  (P),  is  to  find 

X..  and  y,  values  that: 
ij  •'k 

m  n  p 

Minimize  Z  I  a .  ,  x .  .  +  I  b,  y,  (i) 

i=l  j=l  k=l  ^ 


m 


subject  to  Z  X , .  =  1 

i=l  ^ 

j=l , • . • , n 

(ii) 

m  n 

Z  Z  d  ,  X  < 
1=1  j=l  - 

s,  Yi, 
k  k 

k=l, . . . ,p 

(iii) 

X . .  =  0  or 

1  for 

all  j  and  j 

(iv) 

y  =  0  or 

1  for 

all  k 

(v) 

ii 


where  i  ,  j  ,  k  are  Indices  for  designs,  activities,  and  facilities, 
respectively:  has  value  1  if  and  only  if  activity  j  uses  design 

1  ,  and  yj^  has  value  1  if  and  only  if  facility  k  is  used.  A  design 

involves  the  use  of  one  or  more  facilities,  and  ':he  same  design  may  be  used 
by  several  activities. 


Problem  (P)  has  the  objective  of  minimizing  the  sum  of  a..'s  — 

the  variable  costs  due  to  the  assignments  of  activities  to  designs,  and 
bj^'s  —  the  fixed  costs  due  to  the  facilities  used.  Constraints  (ii) 


and  (iv)  ensure  that  each  activity  is  assigned  to  c  single  design, 
d.  is  the  capacity  required  at  facility  k  if  activity  j  uses 

IJK. 


Each 

design 


i  ,  and  is  thus  equal  to  zero  if  design  i  does  not  involve  the  use  of 
facility  k  .  Constraints  (iii),  therefore,  ensure  that  for  each  facility 
k  used,  the  total  capacity  required  does  not  exceed  the  capacity  available 
s^  .  The  difficulty  in  solving  problem  (P)  stems  from  the  indirect  relation¬ 


ship  between  the  assignments  and  facilities,  i.e.,  an  assignment 
bears  on  all  the  constraints  (iii)  for  which  positive, 
therefore,  on  several  y^^  variables. 


"ij 

and , 
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The  branch-and-bound  solution  algorithm  uses  Lagrangian  relaxation 
as  a  basic  step  in  obtaining  lower  bounds.  In  addition,  it  includes  several 
operational  rules,  such  as  a  branching  rule  for  a  judicious  choice  of  the 
branching  variable,  a  capacity  rule  to  eliminate  infeasible  assignments, 
and  a  bounding  rule  to  eliminate  non-optimal  assignments. 


This  dissertation  includes  relevant  background  leading  to  the 
formulation  of  problem  (P) ,  mathematical  development  of  the  branch-and- 
bound  solution  algorithm,  a  detailed  test  example,  and  computational  test 
results  using  the  computer  program.  The  areas  of  application  are  identified, 
and  suggestions  for  further  Improvement  of  the  branch-and-bound  solution 
algorithm  are  included. 

The  computer  program  has  been  written  in  FORTRAN  IV.  A  detailed 
description  of  the  computer  program  and  guidelines  for  its  use  are  Included 
in  a  separate  document  entitled  "Program  Description  and  User's  Guide  for 
ZIPCAP — a  Zero-one  Integer  Program  to  solve  multiactivity  multifacility 
Capacity-constrained  Assignment  Problems."  Although  developed  for  capac¬ 
itated  problems,  the  computer  program  can  also  be  used  to  solve  uncapacitated 
problems  in  which  it  is  assumed  that  the  facilities  have  infinite  capacity. 
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1.  INTRODUCTION 


Muitiact  ivity  mul t ilac il  it y  assignment  p?'obleins  arise  in  such 
diverse  areas  as  public  health  care  systems  and  private  multi-echelon 
inventory/distribution  systems.  Such  systems  involve  the  assignment  of 
activities  or  tasks  to  groups  of  facilities  in  such  a  way  that  total 
system  cost  is  minimized.  The  total  system  cost  has  components  (fixed 
costs)  tliat  depend  on  the  facilities  actually  used  as  well  as  components 
(variable  costs)  that  depend  solely  on  the  assignment  made.  Most  recently 
[Gross,  Pinkus,  and  Soland  (1979)]  there  has  been  interest  in  including 
facility  capacity  constraints  as  well.  For  this  kind  of  problem,  i.e., 
a  multiactivity  multifacility  capacity-constrained  0-1  assignment  problem, 
we  have  developed  a  solution  algorithm  of  the  branch-and-bound  type  and 
a  computer  program  based  upon  it. 

The  computer  program  and  guidelines  for  its  use  are  described  in 
a  separate  document  [Chhabra  and  Soland  (1980)]  titled  "Program  Descrip¬ 
tion  and  User's  Guide  for  ZIPCAP  —  a  Zero-one  Integer  Program  to  solve 
multiactivity  multifacility  Capacity-constrained  issignment  Problems." 

This  document  describes  the  development  of  the  solution  algorithm 
and  computational  test  results  using  the  computer  program.  Suggestions 
for  further  improvement  in  the  solution  algorithm  are  e'.lso  included. 

This  chapter  reviews  tlie  relevant  literature,  provides  background 
leading  to  the  mathematical  formulation  of  the  multiactivity  multifacility 
capacity-constrained  0-1  assignment  problem,  called  problem  (P) ,  and 
includes  potential  areas  of  application.  The  theoretical  base  for  develop¬ 
ing  the  algorithm/methodology  are  described  in  Chapter  2.  Various  components 
of  the  methodology  are  covered  in  detail  in  Chapter  3.  Chapter  4  provides 
an  overview  of  the  computational  procedure  and  the  computer  program,  whereas 
computational  test  results  are  given  in  Chapter  5.  Suggestions  for  further 
research  and  potential  improvements  in  the  algorithm  are  included  in  Chapter 
6. 
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It  may  be  noted  that  the  basic  terminoloi,y ,  described  below, 
in  the  formulation  of  problem  (P)  Includes:  actlv\tles  that  must  be 
assigned,  facilities  which  serve  the  activities,  designs  involving  one 
or  more  facilities,  fixed  costs  associated  with  tne  facilities,  and 
variable  costs  associated  with  the  assignment  of  activities  to  designs. 

The  following  review  of  the  relevant  literature  starts  with  the 
classical  assignment  problem  and  leads  to  the  fonnulation  of  problem  (P) . 
Different  authors  have  used  various  terminologies  in  describing  relevant 
formulations.  In  the  following  discussion,  the  original  terminologies 
are  used,  and  are  followed  by  our  equivalent  terminology,  where  appro¬ 
priate,  shown  in  parenthesis. 

1 . 1  Generalized  Assignment  Problem 

In  a  classical  assignment  problem  [Hillier  and  Lieberman  (1980)], 
the  purpose  is  to  find  optimal  pairs  of  agents  ard  tasks  or  activities. 
Each  task  is  assigned  to  a  single  agent,  and  each  agent  is  given  a 
single  task,  and  the  suitability  of  a  particular  set  of  assignments 
is  determined  by  a  single  criterion  function  such  as  minimization  of 
cost.  In  a  generalized  assignment  problem  (GAP),  several  tasks  can  be 
assigned  each  agent,  subject  to  the  resources  available  to  the  various 
agents  [Ross  and  Soland  (1975)],  e.g.,  assigning  software  development 
tasks  to  programmers  and  assigning  jobs  to  computers  in  a  computer 
network. 

A  variety  of  well-known  facility  location  and  location-alloca¬ 
tion  problems  have  been  shown  to  be  equivalent  to,  and  therefore 
solvable  as  GAP's  [Ross  and  Soland  (1977)].  Here,  in  general,  the 
tasks  represent  demand  centers  for  a  good  or  service,  and  the  agents 
represent  supply  centers  to  be  established  at  potential  sites  or 
locations.  Each  demand  center  must  be  supplied  from  a  supply  center. 

A  fi.xed  cost  is  incurred  for  each  supply  center  established,  and,  in 
addition,  there  is  a  cost  incurred  for  each  unit  processed  at  a  supply 
center  and  transportation  costs  incurred  for  the  units  sent  from  supply 
centers  to  demand  centers.  The  problem  may  be  "uncapacitated"  —  when 
there  is  no  limit  to  the  number  of  units  that  may  be  processed  by 


a  supply  oi  "cap.K- itaied"  —  when  there  are  restrictions  on  the 

number  oi  units  that  may  be  processed.  The  objective  is  to  select  supply 
center  locations  and  set  up  a  distribution  assignment  so  that  the  total 
cost  is  minimized. 

1 . 2  Mult iac tivity  Mul t i t ac il  ity  Uncapacitated 
Assignment  iroblem 

A  salient  feature  of  the  above  facility  location  problems  is  that 
each  demand  center  (activity)  is  assigned  to  a  single  supply  center 
(facility).  Sometimes,  however,  it  may  be  desirable  to  assign  an  activity 
to  more  than  one  facility.  This  leads  to  the  concept  of  "design,"  and 
the  multiactivity  multifacility  assignment  problem  [Pinkus,  Gross,  and 
Soland  (i97J)].  Before  describing  such  a  problem,  some  terminology  is 
considered  first. 

A  design  Involves  the  use  of  one  or  more  facilities,  and  represents 
a  meaningful  configuration  of  facilities  along  with  a  meaningful  strategy 
for  using  them  —  as  Illustrated  in  the  following  examples. 

Consider  five  facilities  and  their  locations  as  shown  in  Figure 
1(a).  (From  practical  considerations,  these  may  be  existing  and/or 
potential  locations.)  Three  of  the  possible  designs  are  shown  in 
Figures  1(b)  to  1(d).  Design  1  is  completely  cent-alized  since  it  uses 
only  one  facility,  whereas  design  3  is  completely  decentralized  since  it 
uses  ail  the  facilities. 


© 

© 

© 

© 

(a) 

Locations  for  five 
facilities 

© 

© 

© 

(a)  Design  2:  three  facili¬ 
ties — locations  ©.  © 
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Figure  1. 
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(b) 

Design  1:  one 

facility- 

locat ion(^ 

© 

© 

© 
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(d)  Design  3:  five  facilities 
all  local  ions 


Examples  of  alternate  designs 
for  a  system  of  five  facilities 
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It  is  [)os;;ible  for  several  designs  to  have  the  same  facilities 
but  different  configuration  and  strategies  for  using  these  facilities, 
e.g.,  a  multiproduct  multi-echelon  inventory  system  [Gross,  Pinkus, 
and  Soland  (1979)].  Figure  2(a)  shows  design  1  containing  certain 
facilities  (warehouses)  at  the  central,  regional,  and  local  levels 
or  echelons.  Figure  2(b)  shows  design  2  with  the  same  facilities 
but  having  a  different  configuration. 
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(a)  Design  1  (b)  Design  2 


Central 

Regional 

Local 


Figure  2. 


Dxap.ple  of  alternative  designs  having 
the  sane  facilities  but  different 
conf  ifiuration 


The  distribution  of  a  given  activity  at  various  facilities  under  design  1 
would  be  different  than  under  design  2,  depending,  of  course,  on  the 
inventory  policies.  This  results  in  different  variable  costs  (described 
later)  for  that  activity  under  design  1  as  against  design  2.  In  fact, 
it  is  possible  to  have  a  situation  where  two  or  more  designs  have  the 
same  facilities  and  the  same  configuration  but  different  strategies, 
resulting  in  different  variable  costs.  For  example,  one  strategy  might 
specify  an  equal  distribution  of  a  specific  activity  over  the  various 
facilities,  whereas  another  strategy  could  impose  a  different  distribu¬ 
tion  scheme  over  the  same  facilities. 

In  general,  if  a  system  is  to  be  composed  of  at  most  p 

facilities,  the  number  of  alternative  designs  is  2*^-1  if  no  two  designs 
have  the  same  facilities.  However,  with  the  same  facilities  but  different 
configurations  and  strategies,  the  number  of  altern-itive  designs  could  be 
much  higher.  In  practice,  it  is  possible  to  eliminate  a  majority  of 
alternative  designs  because  of  geographical,  political,  economical,  and 
other  factors. 


-  5  - 


i'l,L‘  :iuil  t  iai’Liv  ity  mu  1 1  il  ao  1 1  Ity  assignment  problem  seeks 
minimization  oi  some  measure  oi  total  system  ctist  such  as,  total 
expected  cost  over  a  given  time  period  or  total  discounted  cost  over 
the  lifetime  of  the  system.  The  system  cost  will  include  investment 
costs  for  building  or  leasing  tlie  system,  operating  costs  for  operation 
and  maintenance  of  the  system,  and  the  costs  for  providing  necessary 
services.  Bor!\  the  investment  costs  and  the  opetating  costs  have 
fi.'red  as  well  as  variable  components  [Ross  and  Soland  (1980)].  The 
fixed  components  include  those  costs  associated  with  the  facilities 
ot  a  given  design  which  are  independent  of  tlic  activities  served. 

Such  costs  are  called  fixed  costs.  On  the  other  hand,  the  variable 
components  and  the  service  costs  include  tliose  costs  which  are 
completely  dependent  on  the  service  demand  of  the  activities  at 
th.e  various  facilities  in  a  given  design.  Such  costs  are  called 
variable  costs.  By  definition,  both  the  fixed  costs  and  the  variable 
costs  are  relative  terms. 

An  equivalent  formulation  of  the  multiact ivity  multifacility 
assignment  problem  defined  by  Pinkus,  Gross,  anc  Soland  (1973)  is  as 
follows. 


Let  a  .  . 
■LJ 


b 


k 


b 


ik 


variable  cost  of  activity  j  using  design  i 
(i“l,...,mj  j“l,*..,n) 

fixed  cost  of  facility  k  (k=l,...,p) 

1  if  facility  k  is  included  in  design  i  , 


0  otherwise. 


Is  defined  as: 
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The  decision  variable  x. . 

il 


=  1  if  activity  j  uses  desit;n  i  , 
=  0  otherwise. 


Then,  the  uncapacitated  assignment  problem  called  problem 
(PU)  if  to  find  values  that: 


(PU)*C 


Minimize 


subject  to 


where 


m  n 

E  I  a..  X..  +  E  b,  u/  E  b.,  Ex 

"J  k=l  li=l  j=l  "j 


i=l  j=l 


m 


x .  .  =  1 


i=l 
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for  j=l, . , , ,n 


X..  =  0  or  1  for  all  i  and  i 
ij 


u(-)  =  0  if  (•)  <  0  , 
=  1  if  (.)  >  0  . 


(1) 

(2) 

(3) 


The  objective  function  of  this  problem  coisists  of  two  distinct 
parts.  The  first  part  represents  the  total  variabi?  cost,  and  the 
second,  the  total  fixed  cost  of  the  system.  Constraints  (2)  and  (3) 
ensure  that  each  activity  is  assigned  to  a  single  design.  Of  course, 
the  optimal  solution  may  involve  the  use  of  more  chan  one  design. 

Problem  (PU)  is  a  0-1  nonlinear  programming  problem  (because  of 
the  step  function  u  ),  and  a  branch-and-bound  algorithm  using  linear 
underestimates  for  the  nonlinear  part  of  the  objective  function  has  been 
described  in  Pinkus,  Gross,  and  Soland  (1973).  A  heuristic  procedure 
for  this  problem  is  given  by  Khumawala  and  Stinson  (1980)  in  an  unpublished 
paper.  This  procedure  is  an  extension  of  some  earlier  work  [Khumawala 
(1973)]. 

1.3  Adding  Capicity  Constraints  —  Problem  (P) 

A  weakness  of  problem  (PU)  is  that  it  assum::s  unlimited  capacity 
available  at  each  facility  in  terms  of  the  activities  using  a  given 
facility.  In  practice,  a  facility  may  not  have  the  capability  to 
serve  every  activity,  and  may  have  restrictions  as  to  the  total 
capacity  available  to  handle  more  than  one  activity. 
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d..j^  =  capacity  required  at  facility  k  for  activity  j 
when  activity  j  uses  design  i  . 

If  design  i  does  not  include  facility  k  ,  then  d.  =  0  for  all  j 

i.i  k 

Define  the  decision  variable  y,  as: 

■'k 

yj^  =  1  if  facility  k  is  used, 

=  0  otherwise  . 

Then  the  assignment  problem  [Gross,  Pinkus,  and  Soland  (1979)], 

called  problem  (P)  is  to  find  x..  and  v,  values  that: 

✓  ij  ■'k 


Minimize  I  E  a.,  x,  +  E  b,  y, 

IJ  IJ  k^k 


subject  to 


E  X. .  =  1 


j=l . n 


'^ijk  ^ij  -  Vk  k=i,...,p 

^ij’  ^k,  “  ®  ^  i>j.k 


Constraints  (5)  of  problem  (P)  ensure  that  the  capacities  available  at 
the  facilities  are  not  violated.  Problem  (P)  is,  thus,  a  raultiactivity 
muitifacility  capacity-constrained  0-1  assignment  problem,  as  compared 
to  problem  (PU)  which  is  uncapacitated.  In  problem  (P) ,  constraints  (2) 
along  with  the  part  of  constraints  (6)  involving  the  x^^^'s  ensure  that 

each  activity  is  assigned  to  a  single  design.  Of  course,  the  optimal 
solution  may  result  in  the  use  of  more  than  one  design. 

For  an  example  of  five,  facilities  and  three  designs  as  shown 

in  Figures  i(b)  to  1(d),  and  four  activities;  the  matrix  [a^.ib,  Id..,] 

ij  '  k'  ijk 

is  as  shown  in  Figure  3. 

1.3.1  Comparison  with  the  uncapac itated  assignment  problem. 

Comparison  of  the  capacitated  problem  (P)  with  the  uncapacitated 
problem  (PU)  shows  that  the  objective  functions  (1)  and  (A)  are  equiva¬ 
lent  and  constraints  (2)  in  each  are  the  same.  Constraints  (5)  serve 
to  Impose  the  capacity  constraints  and  at  the  same  time,  for  a  given 
design,  the  relevant  facilitie.s  are  forced  in  the  solution.  For  an 
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X  equal  to  1  ,  all  the  facilities  with  d .  . ,  >0  .nust  have  y, 

ij  ^  ijk  ■'k 

values  equal  to  1  in  order  to  satisfy  (5)  and  the  corresponding  fixed 

costs  are  therefore  included  in  (4).  If  ~  ®  *^ijk  ^  ^ 

then  X..  must  be  0  in  order  to  satisfy  (5). 


Problem  (P)  has  been  formulated  as  a  0-1  linear  programming 
problem  wh.ereas  problem  (PU)  was  formulated  as  a  0-1  nonlinear 
programming  problem. 


case  of 
design 


Note  that  problem  (PU)  can  be  easily  obtained  as  a  special 

problem  (P)  by  letting  d...  equal  1  (for  all  j  )  if 

1.)  k 

i  uses  facility  k  ,  and  setting  all  s^^  equal  to  n  . 


In 


other  words,  the  corresponding  formulation  is  to  find 
values  that: 


X  .  . 
11 


and 


y 


k 


(PI) 


m  n 


Minimize  L  L  a..x.,  +  Z  b,  y, 

■  11  tJ  iJ  II  k  K 
1=1  j=l  k=l 


m 

i=l 


subject  to  ^  x^j  =  1 


j=l, . . . ,n 


m  n 

I  e . ,  Z  X . .  <  ny,  k=l , .  — p 
ik  ,  ,  ij  -  -^k 


•  1  A  tv  ,  , 

1=1  j=l 


x.,,y,  =  0  or  1  for  all  i,j,k 
13’ -^k 


^where 


'ik 


1  if  design  i  uses  facilitv  k  , 
0  otherwise  . 


(4) 


(2) 

(7) 

(6) 


1.3.2  Comparison  with  the  fixed-charge  location-allocation  problem. 

Problem  (P)  bears  a  resemblance  to  the  well-known  fixed-charge 
location-allocation  problem  or  capacitated  facility  location  problem 
(Ceoffrion  (1975);  Ross  and  Soland  (1977)].  Tliere  are,  however, 
very  significant  differences  between  the  two.  In  order  to  point  out 
these  differences,  here  is  a  statement  of  the  location-allocation 
problem  (Ij\)  as  given  by  Cross,  Pinkus,  and  Soland  (1979)  in  a  form 


_  10 


similar  to  that  oi  problem  (P) . 


Find  and  yj^  values  tliat 


(LA)  -C 


Minimize 


subject  to 


PR  P 

ii  j'l  ‘“j  ’‘“j  ^ 


I.  X  =  i 

k-1  '‘J 


Jl  “j  "kj  -  ^''k 


j“i» • • • >R> 


k=l , . , 


>  0  ,  yj^  =  0  or  1  Jor  all  j 
.ind  k 


(8) 


(9) 


(10) 


(11) 


Here  represents  the  fraction  of  customer  (acti^■ity)  j's  demand 

that  is  supplied  by  a  facility  at  location  k. 

The  most  important  distinction  between  problem  (LA)  and  problem 
(P)  is  the  relationship  between  assignments  and  facilities.  In  problem 
(LA)  there  is  a  direct  connection  between  the  assignments  made  and  the 
facilities  required,  and  each  assignment  affects  cnly  one  facility,  i.e., 
the  assignment  >  0  has  a  bearing  on  only  one  of  the  constraints 


(10)  and,  therefore,  on  only  one  variable  y 


On  the  other  hand , 


in  problem  (?),  the  connection  between  the  assignments  made  and  the 
facilities  required  is  indirect,  and  each  assignment  can  affect  several 
facilities,  i  e.,  the  assignment  x^^  =  1  bears  on  all  of  the  constraints 

(5)  for  which  d.  >  0  and,  therefore,  on  several  variables  y,  . 

ij  K  K, 

Another  distinction  is  the  relative  difficulty  of  the  two  problems. 
While  problem  (LA)  is  not  easy  to  solve,  brancl»-and-bound  approaches 
have  been  successful  in  dealing  with  it  because  oncz  values  are  specified 
for  the  y^^  ,  the  x.^^  are  found  by  solving  a  transportation  problem. 

Problem  (LA)  becomes  more  difficult  if  the  constraints  x  >  0  in  (10) 

kj  - 

arc  replaced  by  x  =  0  or  1  in  order  to  preclude  supply  of  customer 
(activity)  j's  demand  by  more  than  one  facility.  With  this  change. 


IJ 


problem  (LA)  may  be  treated  as  a  generalized  assignment  problem  and 
is  solvable  using  an  efficient  branch-and-bound  algorithm  [Ross  and 
Soland  (1977)].  Problem  (P)  is  more  difficult  than  this  variation 
of  problem  (LA)  because  of  the  above  stated  indirect  connection 
between  the  assignments  and  the  facilities.  Even  after  values  have 
been  specified  for  ail  the  ,  problem  (P)  remains  a  difficult 

0-1  linear  programming  problem  because  of  the  interaction  of  the 
constraints. 

1.3.3  Solving  problem  (P). 

The  capacitated  problem  (P)  has  mn+p  0-1  variables  and 

n+p  constraints,  so  the  problem  dimensions  may  be  large  from  practical 
considerations.  For  example,  with  m=n=30  and  p=20  ,  problem  (P) 
has  920  variables  and  50  constraints.  The  0-1  LP  computer  codes 
generally  available  are  limited  in  terms  of  problem  size.  For  example, 
the  code  used  by  Gross,  Plnkus,  and  Soland  (1979)  can  handle  up  to  AO 
variables  and  20  constraints.  A  better  and  more  efficient  code 
[Geoffrion  and  Nelson  (1968)]  allows  up  to  90  vaviablv=!S  and  50 
constraints.  This  fact,  together  with  the  structure  of  problem  (P) 
suggests  that  a  specialized  algorithm  could  be  developed  that  would 
be  more  efficient  for  practical  problems  than  the  general  integer 
linear  programming  algorithms  (on  which  the  available  codes  are 
based) . 

With  the  above  background  in  mind,  the  development  of  the 
solution  algoritlwi  and  the  computer  program  to  solve  problem  (P) 
was  undertaken  and  is  described  in  Chapters  2  through  A. 

1 . A  Areas  of  Appl ication 

The  solution  algorithm  and  the  computer  program  are  designed 
to  solve  a  multiactivity  mviltifacility  capacity-constrained  0-1 
assignment  problem,  i.e.,  one  which  can  be  formulated  as  problem  (P) . 

The  basic  elements  of  such  a  problem  are  activities  that  must 
be  assigned,  facilities  and  their  meaningful  configurations 


represented  as  designs,  the  fixed  and  variable  costs,  and  the 
capacity  requirements  of  the  activities. 


The  formulation  (P)  applies  to  existing  a.ad/or  proposed 
facilities.  In  other  words,  it  is  useful  for  a  situation  where  the 
decision  may  be  to  delete  some  of  the  existing  facilities,  as  well  as 
for  a  situation  where  the  decision  involves  a  selection  out  of  a  set 
of  proposed  facilities. 

Table  1  includes  examples  of  areas  where  formulation  (P)  is 
applicable.  Witliin  each  application  area,  activities  and  facilities 
are  identified.  The  implications  of  designs,  variable  costs,  and 
fixed  costs  are  apparent. 

Obtaining  the  values  of  the  data  elements  b,  ,  d .  ,  ,  s  , 

K  ijk  k 

and  in  particular  a_  ,  can  be  a  simple  or  a  complex  exercise  depending 

on  tl'.e  particulars  of  the  application,  and  the  nature  of  the  components 
comprising  these  elements.  For  example,  in  designing  multi-echelon 
inventory  sy3teus  [Gross,  Pinkus,  and  Soland  (1979)],  a^^  represents 

the  inventory  cost  of  product  (activity)  j  using  echelon  structure 
(design)  i  and  bj^  represents  the  fixed  cost  C’f  installation  (facility) 

k  .  The  inventory  cost  a^^  includes  the  cost  of  procurement,  carrying 

inventory,  filling  orders,  and  stockouts.  The  value  a^^  ,  and  associated 

inventory  stockage  policies,  are  arrived  at  by  solving  a  multi-echelon 
inventory  problem.  In  other  words,  for  product  j  stocked  under  echelon 


structure  i 
Installat ion 


,  optimal  inventory  policies  are  determined,  at  each 
of  the  structure,  which  yields  a^^  .  The  facility  fixed 


cost  bj^  includes  the  capital  expenditure  for  building  the  installation. 


along  with  a  number  of  fixed  costs  associated  with  operating  it,  such  as 
administrative  expenses,  the  expense  of  renting  the  facility  (if  it  is 
not  built),  .nd  certain  other  fixed  operating  expenses. 


In  the  case  of  designing  a  support  system  for  repairable 


EX.MU'LKS  of  application  areas 


*Gross  and  Pinkus  (1979) 
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items  [Gross  and  Pinkus  (1979)],  represents  the  total  variable 

cost  if  unit  type  (activity)  j  is  repaired  under  design  i  .  The 
set  of  parameters  taken  into  consideration  to  compute  this  cost  for 
each  unit  type  includes  such  things  as  varying  population  sizes, 
failure  rates,  average  repair  times,  costs  associated  with  their 
repair,  the  purchase  and  storage  of  spares,  the  purchase  of  repair 
channels,  and  travel  to  depots  (facilities)  for  repair.  A  computer 
program  is  used  to  solve  a  spares  and  server  provisioning  problem,  and 
the  results  provide  the  basic  information  to  compute  a^^  . 

Thus,  in  general,  the  data  elements  of  problem  (P)  may  be  obtained 
directly  and/or  by  solving  other  related  problem(s);  it  depends  on  the 
definition  and  the  nature  of  the  components  comprising  these  data  elements 
for  a  specific  application  area. 


2.  DEVELOPMENT  OF  THE  SOLUTION  ALGORITIIM 


The  solution  algorithr'.  that  has  been  developed  ti  solve  problem  (P) 
is  a  branch-and-bound  procedure  which  makes  use  of  Lagrangian  relaxation 
as  a  basic  step. 

This  chapter  considers  two  different  Lagrangian  relaxations  of 
problem  (P),  their  general  characteristics,  and  some  useful  results  leading 
to  the  specific  case  of  Lagrangian  relaxation  utilized  in  the  solution 
algor  itlim. 

2 . 1  Lagrangian  Relaxation 

Taking  a  set  of  "complicating"  constraints  of  a  ganeral  mixed- 
integer  program  into  the  objective  function  in  a  Lagrangian  fashion  (with 
fixed  multipliers)  results  in  a  "Lagrangian  relaxation"  of  the  original 
problem  [Geoff rion  (1974)].  The  relaxed  problem  is  easy  to  solve  compared 
to  the  original  problem,  and  provides  a  lower  bound  (for  minimization 
problems)  on  the  optimal  value  of  the  original  prob-.em. 

Although  the  use  of  Lagrangian  relaxation  in  discrete  optimization 
has  been  reported  prior  to  1970  [e.g.,  Lorie  and  Savage  (1955),  Everett 
(1963),  and  Gilmore  and  Gomory  (1963)],  the  "birth"  of  the  Lagrangian 
approach  as  it  exists  today  [Fisher  (1978)]  occurred  in  1970  with  the 
successful  application  of  Lagrangian  relaxations  to  the  traveling  salesman 
problem  [Held  and  Karp  (1970,  1971)],  This  was  followed  by  application 
or'  Lagrangian  relaxation  to  scheduling  problems  [Fisher  and  Schrage  (197  2),  and 
Fisher  (1973,  1976)],  the  general  integer  programming  problem  [Shapiro 
(1971),  and  Fisher  and  Shapiro  (1974)]  and  the  generalized  assignment  problem 
[Ross  and  Soland  (1975)].  Table  2  lists  the  applications  of  Lagrangian 
relaxation  as  given  by  Fisher  (1978).  A  review  of  Lagrangian  relaxation 
is  also  provided  by  Shapiro  (1977)  and  Christofidcs  (1930). 
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TABLE  2 

* 

APPLICATIONS  OF  LAGRANGIAN  RELAXATION 


Problem 

Researchers 

Lagrangjan  Problem 

TRAVELING  SALESMAIJ 

Symmetric 

Asymmetric 

Symmetric 

Asymmetric 

SCHEDULING 

n|m  Weighted 
Tardiness 

1  Machine  Weight 
Tardiness 

Power  Generation 
Systems 

GENER.\L  IP 
Unbounded  Variables 
Unbounded  Variables 
0-1  Variables 

LOCATION 

Uncapacitated 

Capacitated 

Databases  in 

Computer  Networks 

GENERALIZED  ASSIGNMENT 


Held  &  Karp  (1970,  1971) 
Bazarra  &  Goode  (1977) 

Balas  &  Christofides  (1976) 
Balas  &  Christofides  (1976) 


Fisher  (1973) 

Fisher  (1976) 

Muckstadt  &  Koenig  (1977) 

Fisher  &  Shapiro  (1974) 
Burdet  &  Johnson  (1976) 
Etcheberry,  et.  al.  (1978) 

Cornuejols,  Fisher,  & 
Nemhauser  (1977) 

Geoffrion  &  McBride  (1977) 
Fisher  6  Kochbaum  (1978) 

Ross  &  Soland  (1975) 
Chalmet  &  Gelders  (1976) 


SET  COVERING— PARTITIONING 

Covering  Etcheberry  (1977) 

Partitioning  Nemhauser  &  Weber  (1978) 


Spanning  Tree 
Spanning  Tree 
Perfect  2-Matching 
Assignment 


Pseudo-Polynomial 
Dynamic  Programming 
Pseudo-Polynomial  DP 
Pseudo-Polynomial  DP 


Group  Problem 
Group  Problem 
0  -  1  GUB 

0-1  VUB 
0  -  1  VUB 

0  -  1  VUB 


Knapsack 

Knapsack,  0-1  GUB 

0  -  1  GUB 
Matching 


*Sourcc:  Fisher  (1978) 


-  17 


2.1.1  Kel.ixuif,  Problem  (P) 

Uv  Jividine,  cuastraints  (5)  by  s  and  letting  ~  •  r  / c  » 

K. 

problem  (P)  can  be  restated  as  follov;s. 


Minimize  Z  Z  a.,  x..  +  Z  b,  y, 

i-1  j=i  ^ 


sub i ec  t  to  I  x  .  .  =  1 
i=i 


j=_, . . . ,n 


^  ^iik  ""ii  -  . . 

i=l  j=l 


-x..,  y,  =  0  or  1  for  all  i,i,k 
ij  k 


A  Lagraagian  relaxation  (LR^)  of  problem  (P)  relative  to  constraints 


(2)  is  obtained  as 


.Minimize  Z  I  a.,  x..  + 

i=l  j=l  k 


(LR^)<  subject  to  _Z^  _Z^r.  x..  <  yj^ 


P  n  /  m  \ 


X .  . ,  y,  =  0  or  1 
ij  ^k 


k=l, .  . . ,p 


for  all  i,j,k 


wliere  the  are  Lagrange  multipliers;  it  follo'js  that  the 

uptim.jl  value  of  (iroblera  (LR^^)  is  a  lower  bound  on  the  optimal  value  of 

problem  (P).  i.e.,  Z(LR  )  <  Z(P)  .  We  will  continue  to  use  this  notation 

u  — 


in  wliich  /.(•)  is  the  optimal  value  of  problem  (•)• 


Another  Lagrangian  relaxation  (LR^)  of  problem  (P) ,  relative  to 


constraints  (5'),  is  obtained  as 
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Minin'  ize 


E  Z  a  X .  ,  +  b  y  -  E 
1  j  «  k  “  k 


"k  (/k  -  ^  5  'ijk) 


subject  to  (2)  and  (6),  or  equivalently. 


(LR  ) 


Minimize 


E  E  X  .  / 

ij 


a. .  +  E  V,  r . . 
ij  k  k  13 


ik)  -  ^  fk  -  \) 


subject  to  E  X..  =  1 


i  "J 


j=l , . . . ,P 


X..,  y,  =  0  or  1  for  all  i,j,k 
ij  R 


(13) 

(2) 

(6) 


where  the  are  non-negative  Lagrange  multipliers;  it  follows 

that  Z(LR  )  <  Z(P)  , 

V  — 


2.1.2  General  Characteristics 

A  Lagrangian  relaxation  provides  a  lower  bound  on  the  optimal  value 
of  the  orginal  problem,  i.e.,  in  our  case  Z(LR^)  ^  Z(P)  and  Z(LR^)  :<  Z(P) 
The  usefulness  of  a  Lagrangian  relaxation  depends  on  the  closeness  of  this 
lower  bound  to  the  optimal  value  of  the  orginal  problem.  However,  the 
relaxation  must  be  "easy"  to  solve  relative  to  the  original  problem.  We 
observe  that  the  optimal  value  of  yj^  in  problem  Is  1  if 

(v,  -  b, )  >  0  and  0  if  (v,  -  b,  )  <  0  ,  and  then  prol lem  (LR  )  reduces 
kk  —  kk—  v 

to  n  0-1  "multiple  choice"  problems  which  are  very  easy  to  solve. 

On  the  other  hand,  problem  (LR  )  reduces  to  k  0-1  knapsack  problems. 

u 

However,  these  problems  are  not  independent  because  of  the  interaction 
of  constraints  (5')  and  the  indirect  relationship  described  earlier  in 
Section  1.3  between  the  assignments  and  the  facilities.  In  view  of  this 
complexity,  relaxation  (LR^)  will  not  be  considered  further. 

The  choice  of  Lagrange  multipliers  in  relaxation  (LR^)  should 

be  such  that  Z(LR^)  is  as  large  as  possible  and  hence  as  close  as 

possible  to  Z(P)  in  view  of  the  relationship  Z(LR^)  ^  Z(P)  .  In  other 

words,  an  equivalent  problem  is  to  find  a  vector  v  ! representing  v^^,  V2t 

to 
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(») 


Miixlmize 

V  ■>  0 


(14) 


Obviously,  X,(1,R^)  <  Z(D)  <  Z(P) 


The  general  properties  ot  Lagrangian  relaxation  have  been  well 
described  in  the  literature  [e.g.,  Ceoffrion  (1974),  Geoffrion  and  McBride 
(197b),  and  Fisher  (1980)].  Some  oi'  these  properties  relating  the 
Lagrangian  relaxation  and  the  usual  LP  relaxation  are  stated  below. 

The  LP  relaxation  (P)  of  problem  (P)  is  obtained  by  relaxing  the 


integrality  constraints 

(6),  i.e. 

,  tlie  formulatioii  (P)  is 

r 

.Minimize 

m 

); 

i=l 

n 

la  X 

j=i 

P 

.  .  +  Z  by 

IJ  k 

(4) 

m 

subject  to 

L'  X  , 
i=l 

1  j=l,...,n 

(2) 

(-”) 

ra 

n 

Z 

i=i 

Z  r. 
j-i 

^ij  -  ^=1 . . 

(5’) 

k=l, .  . . ,p 

(15) 

"ij’ 

>'i.i  ° 

for  all  i , j , k 

(16) 

Note  that  the 

constra ints 

x^^  <  1  are  implicit  in  constraints  (2) 

Also  consider 

the 

following  partial  convex  hull  relaxation 

(P*) 

of  problem  (P) . 

< 

m 

n 

P 

Mini mize 

Z 

i=l 

Z  a  .  .  X 

.1=1 

.  .  +  Z  by 

u  k-1  '' 

(4) 

(!’*)  . 

subject  to 

m 

Z 

i=l 

n 

Z  r  ,  . ,  X 
j=i 

U  '’k  . . 

(5’) 

X  .  . 

ij 

,  y^^  £  convex  iuill  {(2),  (6)} 

(17) 
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Tlitin  till'  relationships  between  the  optiiial  values  of  various 
problems  [Geoff rion  and  McBride  (197B)]  are  as  follows. 

Zip)  j-  Z(1.R'^)  <  max  Z(LR  )  =  Z(D)  Z(:>*)  <  Z(P)  (18) 

V  >  0 

%  a.  %  'V 

where  V  are  the  values  v^,  V2,-..,Vp  of  a  dual  optimal  solution 
of  (P)  corresponding  to  constraints  (5'). 

Tl'.us,  the  optimal  dual  solution  associated  with  the  usual  LP  relax- 
atior.  furnishes  a  choice  of  Lagrange  multipliers  sucii  that  the  associated 
Lagrangian  relaxation  is  at  least  as  tight  as  the  usual  LP  relaxation, 
and  generally  a  good  deal  tighter  and  even  as  tight  as  the  partial  convex 
hull  relaxation. 

Since  Z(D)  =  Z(P*)  ,  the  quality  of  the  bound  obtained  from  the 

Lagrangian  relaxation  depends  on  where  Z(P*)  lies  in  the  range  between 

Z(P)  and  Z(P)  .  It  turns  out  that  problem  (LR  )  possesses  the 

\ 

"integrality  property,"  i.e.,  the  optimal  value  of  problem  (LR^)  is  not 

altered  by  dropping  the  integrality  conditions  on  its  variables  and 
therefore  [Geoffrion  (197A)] 

Z(D)  =  Z(P*)  =  Z(P)  (19) 

ilius,  the  Lagrangian  relaxation  (LR^)  is  no  better  than  the  LP 

relaxation  (P)  .  On  the  other  hand,  Lagrangian  relaxation  does 

not  possess  the  integrality  property  and,  hence,  couid  provide  an  equal 
or  better  bound  tiun  the  LP  relaxation  (P)  ;  but  the  computational 
difficulties  do  not  favor  pursuing  formulation  • 

It  is  possible  to  consider  alternative  formulations  of  problem 
(P)  with,  the  objective  of  obtaining  tighter  bounds.  This  aspect  is 
discussed  in  Ch.ipter  6. 
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2.2  Sumo 

iVo  now  turn  to  the  basic  question  of  choosing  Lagrange  multipliers 
V  so  that  (LK  )  is  optimal  to  the  extent  possil  lo,  \.’hicii  is  equivalent 
to  solving  problen!  (D)  .  We  also  need  to  consider  thrs  questioi’  when  !;ome 
of  the  X..  and  y,  variables  have  been  assigned  values  of  1  or  0,  i.e., 

IJ  R 

at  a  node  other  ti'.an  the  starting  or  "root"  node  in  the  branch-and-bound 
tree.  For  this  purpose,  some  terminology  is  defined  and  formulations 
corresponding  to  problems  (P),  (LK  )  and  (D)  are  first  developed.  Then 

some  important  results  pertaining  to  the  choice  of  Lagrange  multipliers 
will  be  proved.  Gavish  (1978)  provides  a  method  of  obtaining  the  'best' 
multipliers,  based  on  solving  an  equivalent  linear  programming  problem. 
Such  a  formulation  is  difficult  in  our  case,  and,  besides,  we  propose 
to  avoid  solving  LP  problems  in  our  branch-and-bound  procedure. 

Define  the  sets 


S  =  {(i,j)jx^j  has  an  assigned  value  of  1  or  0}  ,  and 
T  =  assigned  value  of  1  or  0)  . 


These  sets  represent  the  partial  solution  of  problem  (?)  and  the  variables 
contained  in  tliese  sets  are  termed  fixed  variables.  [Geoff rion  (1967)]. 

Let  S  and  T  represent  the  corresponding  complementary  sets,  i.e.. 


comprised  of  t!ie  x. 


and 


yj^  variables,  which  have  not  been  assigned 


specific  value.s  .nd,  tiierefore.  are  called  free  variables.  A  completion 

of  a  partial  solution  is  defined  as  a  solution  that  is  determined  by 

S  and  T  together  witii  a  binary  specification  (0  or  1)  of  the  values 

ol  the  free  x..  and  y,  variables  from  sets  S  and  T  , 

11 


Let  SUS  =  Sj  and  TUT  =  Tj  . 

Consider  a  partial  solution  to  problem  (P)  in  which  specific  values 
(of  1  or  0)  are  .issigiied  to  some  of  the 


X..  and  y,  such  that 
ij 
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m 

E  X  <  1  Vj  , 

i=l  ^ 

(i.j)eS 

and  Z  Z  r.  x..  <  y,  VkeT 

Ilk  11  —  k 

1  j 

(l,j)eS 

Z  Z  r .  X .  .  <  1  VkeT 

i  j  - 

(i.j)eS 


ana  sue 


.;h  that  x..  =  1  and  e  =  1  imply  tiiat  keT  and  y,  =  1. 


Recall  that,  by  definition,  -  1  if  design  i  uses  facility  k  , 

and  e  ,  =  0  otherwise, 
ik 

The  problem  of  finding  an  optimal  completion  of  the  partial 
solution  of  problem  (P)  can  be  stated  as  follows. 


Minimize 


Z  Z  a,,  X,,  +  Z  b.y.  +  Z  Z  a^,  x, ,  +  Z  b.y 


1  J  _ 
(i,j)GS 


ij  ij 


I  j  ij  ij  k  (20) 


keT  (i,j)eS 


subject  to 


Z  X,  .  =  1  -  Z  X 
1  ij  1  ij 

(i,J)es 


^  J  'ijk  *1J  i  >'k  -  “  'ijk  ’’fj 


(i,j)eS 


(i,j)eS 


X . . ,  y,  =  0  or  1  V(i,i)eS  ,  kef 

Ij  k 


We  call  this  problem  (Pj^)  where  i,  indicates  the  nodr  in  the  branch-and- 
bound  tree. 


A  Lagranglan  relaxation  of  problem  (P^^)  with  respect  to  constraint: 
(22)  is  obtained  by  introducing  non-negative  Lagrange  multipliers  Vj^  , 


k*l,2,...,p  ;  the  relaxation  is  then 
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Minimi zn 


-  a  .  .  X .  .  +  b  y  +  L  X  a ,  .  x .  .  +  Z  b,  y, 

i  j  _  ij  ij  k  k  .  .  ij  ij  k^'k 


(i,j)eS 


keT 


Zv  y  -  ZZ  r.  x,.-ZZ  r.  x.. 

,  k  k  .  .  ijk  ij  .  .  i]k  ij 

(l,j)es  (l,j)eS 


sub  i  ec  t  to  Z  X  .  .  =  1  -  Z  x . 

i  ""  i  ' 

(i,j)£;S  (i,j)eS 


X .  . ,  y  =  0  or  1 
ij  •'k 


V(i,i)tS  ,  keT 


Rearranging  (24),  and  using  the  reJationship  =  TUT,  we  have  problem 


'■■S.v*  ^ 


1  Minimize 


(LR  ) 

X,,  V 


Z  Z  X  /a  .  +  Z  V.  r .  \  +  Z  Z  X  .  .  I  a .  ,  +  Z  v,  r  .  | 

IJ  (  ‘J  kcT,''  ‘J")  1  j  Ul  U  k  ,3^ 

(i,j)£S  \  ^  '  (ij)es  '  ^ 

'keT  ^  'k^T 


subject  to  Z  X..  =  1  -  Z  X.. 


(1, j)eK 


X .  ,  ,  y  =  0  or  1 
ij  -^k 


V(i,j)eS  ,  keT 


Then  we  iiave  Z(LR,^  ^  important  probieni  is  the  choice  of 


Lagrange  multipliers  v^,v^,...,  ,  represented  by  vector  v  ,  that 

maximize  Z(LRj^  ,  l.e.,  the  problem  (Dj^)  : 


Maximize  |^Z(LK^ 


V  '•  0 
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We  now  state  and  prove  some  theorems  related  to  the  choice  of  Lagrange 
multipliers  .  . .  .v^  . 


Theorem  1:  There  exists  an  optimal  solut  .on  to  problem  (D)  in 
which  ^  for  all  k  . 


Proof : 


Suppose  ^  ^1  *  optimal  solution  to  problem 

(D),  i.e,,  Z(D)  =  ZCLR  .)  where  /  *  <  b,  . 

1  X 


Recall  that 

Z 


(LR  *)  =  Min  Z  Z  X ,  .  f  a .  ,  +  Z  V*  r .  \  y,  /v*  -  b,  ^ 

V*  i  .  11  11  k  ijk^  V 


s.  t .  Z  X ,  .  =1 
i 


Vi 


^ij’  ~  ^  ^  Vi.j.k 


(2) 

(6) 


For  <  bj^  ,  the  optimal  value  of  is  0  ,  and  the  term  -  y^^  (v*  -  bj^ 

in  the  objective  function  is  0. 


Consider  what  happens  if  we  increase  v*  to  .  Call  the 

resulting  vector  ^  .  Consider  problem  (LR^)  •  The  optimal  value  of 

y^^  in  problem  (LR^)  is  0  or  1  ,  and  the  term  ~  *— 1  ~  '’l^  ^  ' 

However,  the  optimal  value  of  y^^  is  the  same  in  pioblems  (LR^*)  and 

(LR  )  for  all  k  >  1  .  Therefore,  the  quantity  Z  y^^  (Vj^  -  b^^)  is  the 

—  k 

same  at  the  optimal  solution  for  both  v  =  v*  and  v  =  v  . 


Since  ^  v*  ,  we  note  that  in  the  objective  function. 


P  P 

a  +  Z  V  r  >  a.,  +  Z  v*  r 

k=l  ~  ^ 


Vi, j  , 


k=l 


and  therefore  Z(LR^)  ^  Z(LR^*) 
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It  tollows  that  thore  is  an  optimal  solut-lon  to  problem  (D) 
in  which  ^  • 

Since  the  choice  of  k=l  was  arbitrary,  the  same  result  holds  for 

any  value  of  k,  k=l,...,p  ;  hence,  there  exists  an  optimal  solution  to 

problem  (D)  in  which  v,  >  b  for  all  k  . 

k  —  k 

Tliocir em _ 'll  There  exists  an  optimal  solution  to  problem  (D^)  in 

v;hich  >  bj^  if  (i)  kcT  or  (ii)  keT  and  y^^  =  0 


Proof : 


Suppose  <  b^  in  an  optimal  solution  to  problem 


(D^)  ,  i.e.,  Z(D^)  =  Z(LRj^  ^*)  where  v*  <  b^ 

Then  k=l  can  he  such  that  kt.T  or  keT  . 


Case  (i):  Let  keT  . 
Recall  that 


Z(LR.  *)  =  Min  11  x..(a,+Ivfr.., 

^.v  i  j  _  \  k 

(i, j)cS 


+  Z  Z  X ,  .  (i 

ij  V 


a.  .  +  Z  V*  r .  . 


13 


(i,J)eS 


“  \ 

keT 


'^k  -  \ 


k  ijk, 


■  ^ 
kex 


s.t.  Zx..=l-Zx.. 

i  ij  i  iJ 

(i,j)cS  (i,j)cS 

X .  . ,  y,  =  0  or  1 
ij’  ^k 


Vi  (21) 

V(i,j)cS  .  kef  (23) 


For  V*  <  bj^  ,  and  kcT  ,  the  optimal  value  of  y^^ 
term  -  y^^  (v*  -  bj^)  in  the  objective  function  is  0  . 


is  0  and  the 


Zb 


Lot  V*  1)0  increased  to  b  ;  call  the  resuxtlng  vector  v 


Consider  probloin  (LR.  )  .  The  optimal  value  of  y,  in  (LR.  )  is 

£,v  1  Ji,v 

0  or  1  ,  then  the  term  -  (v^  -  b^)  is  0  .  For  k  >  1  ,  the 


optimal  value  of  y  being  the  same  in  (LR,  *)  and  (LR.  )  ,  we 
k  x.,v  Jl,v 


find  that  Z  y  (v  -  b  )  is  the  same  at  the  oi timal  solution  for  both 


V  =  V*  and  v  =  .  But  ^  v*  ;  therefore 


^k"ijk^"ij  V.  ''k'ijk  V(i,j)es  and  (i.j)£S 

A.  ^ 


Hence,  Z(LR  )  >  Z(LR.  *)  ,  wherefrom  it  follows  that  there  is  an 
x.,_v  —  i,,v 


optiraal  solution  to  (D.)  in  which  v  >  b  .  S'.nce  the  choice  of 


k=l  was  arbitrary,  the  same  results  hold  for  any  value  of  k,  kef  . 


Hence,  there  exists  an  optimal  solution  to  problem  (Djj )  in  which 


v.  >  b,  for  all  keT  . 
K  —  k 


Case  (ii):  Let  keT  and  =  0 


Considering  problem  (LR,  ■;<)  ,  for  k  =  1  ,  v*  <  b  ,  and 

^  j  V  X  X 


=  0  ,  th.e  term  -  Yj^Cv*  -  b^)  in  the  objective  function  is  0  . 

Increase  v^^  to  b^^  and  call  the  resulting  vector  y_  .  The 
term  -  Y^^CXj;  ■"  ^  1  >  the  optimal  values  of  y^^ 


are  the  same  in  problems  (LR,  *)  and  (LRp  )  .  Therefore  Z 

X.,v  )t,v 


Yj^CVj^  -  b.^)  is  the  same  at  the  optimal  solution  for  both  v  =  v* 


1 


and  V  =  y  .  Since  ^  ’ 
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Therefore  Z(LR,  )  >  Z(LR.  *)  .  It  follows  that  there  exists  an  optimal 

i,v  -  I, V 

solution  to  iti  which  choice  of  k=l  being  arbitrary, 

tlie  same  results  hold  for  any  value  of  k  ,  keT  and  ’  which  proves 

case  (ii)  of  the  Tlieorem. 

It  may  be  added  tliat  there  is  anotlier  possibility  which  complements 
case  (ii)  of  Theorem  2,  i.e.,  if  kcT  and  y^^  =  1  .  We  treat  this  possi¬ 
bility  as  a  conjecture  since  a  result  similar  to  the  one  above  could  not 
be  proved,  as  discussed  now. 


With  ^  '•’f  ■■  ''2^  »  '‘'•2  observe  from  problem  (LK,.  ^*) 

that  for  a  solution  vector  X”  (witli  elements 
elements  y^* ,  .  •  .  ,  y'"'' !  y^*"  1  y2*»‘**>y  *  =  0  or  1)  , 


Z(LR,  *)  =  EEx,.*  /  a. .  +  V*  r. ..  +  Z  v*r..,\ 

^.v  i  j  V  1  k>l 

(i, j)cS 

+  ZZx,,*  (a..+v*r,  +Z  v,  *r.,,\ 
i  j  ij  \  ij  1  iji  i,>l  ^ 

(i,j)eS 


Since  v^*  <  b^  and  y^^  =  1  the  term  -  ”  ^1^  positive. 

If  we  raise  v^*  to  b^  ,  say  v^^  ,  the  term  -  y^^X^  ~  ®  • 

The  difference  between  Z(LR,  ^*)  and  the  objective  function 


problem  (LR^ 

,V> 

with 

X=X  and 

Y=Y  is 

Z  Z  X.,*  V* 

+  (b, 

-  v/')  -  Z 

Z  X  *  b 

i  j  ' 

Ijl 

1 

1 

1 

j  'j 

(b^  -  v^v..)  _ 

Z  Z 

i  j 

X  .  .  * 

iJ 

-  ''l*) 

’’iji  . 

This  difference  can  be  eitlier  negative  or  positive,  and  s )  we  cannot 
conclude  that  there  is  an  optimal  solution  to  problem  (Dj,)  in  which 


.  ^  ^  <g*.  « •- W  »-  ^ 
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.  Wo  boiiove  this  CDiiciusion  to  be  false. 

ic  * 

Theorem  3:  Lot  (X  ,  Y  )  solve  problem  (LR  )  for  v,  =  b, 

- — "" — ^ —  V  k  k. 

*  * 

for  all  k  .  If  (X  ,  Y  )  is  feasible  for  problem 

(P) ,  there  exists  an  optimal  solution  to  problem  (D) 

in  which  v,  =  b,  for  all  k  . 
k  k 

Proof :  In  view  of  Ttieorem  1,  there  exists  an  optimal  solution 

to  (D)  in  which.  Vj,  ^  bj^  for  all  k  ,  i.e.,  v  >  b  . 

Let  ;v  be  such  an  optimal  v  .  We  will  show  that 
Z(LR^)  ^  Z(LR^)  ,  from  which  it  follows  that  v  =  b 

solves  problem  (D)  . 


Recall  that 


Z(LR  )  =  Min  Z  Z  x. . 


X,Y  i  j 


ij 


s.  t,  Z  X  .  .  =1 


ij 


X . 
ij 


yk 


=  0  or  1 


Vj 

Vi,  j,k 


(2) 

(6) 


Since  v>b,  y  =  lVk  is  an  optimal  choice. 
-  k 


Hence,  Z(LR  )  =  liin  Z 

V 

-  X  i 


Z  Z  X .  .  I  a .  .  +  Z  v  r  I  -  Z!  /  v  -  b,  J 
i  j  ^  ij  ijnj  kj 


s,  t , 

Z  X  .  =  1 
i 

Vj 

(2) 

X  =  0  or  1 

ij 

Vi,j 

(6a) 

Now  consider  - 

Since  V  =  b  , 

the  last  term  of  the  objective 

function  drops  out,  and  we  have 


Z(LR,  )  =  Mi 
^  X 


in  Z  Z  X  I  a  +  Z  b,  r ,  j 

,Y  i  j  iJ  \  iJ  k  "J7 


subject  to  (2)  and  (6) 
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-  Min  ;;  X.. 

( a  ,+  y.  b  r 

X  i  j 

y  ij  k  ijk 

subject 

to  (2)  and  (6a) 

Z  Z  X  .  .  *  ^ 

+  Z  b  r  J 

ii 

ij  k  ij*"/ 

where  x’*  with,  elements  x..*  is  the  mirimizinc  solution  vector 

ij 

which  satisfies  (2)  and  (6a)  . 


"k  k 

Now  (X  ,  Y  )  feasible  for  (P)  implies  that 


subject  to  (2)  and  (6a) 


subject  to  (2)  and  (6a) 


_  ZZx..*  ( a . .  +  Z  r .  b,  \ 

1  j  iJ  ^  iJ  k 

,  (-k  K)ij  O  Ok  ,  (-k  k) 

<  ZZx..*  /a..+Zr..,  b  \  =  Z(LR,) 

-  i  j  11  ^  iJ  k  V  ^ 


.  ■A:S-  ■  -’-k  . 
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by  (27),  or  Z(LR^)  ^  Z(LRj^)  ;  it  follows  Chat  v  =  b  solves 
problem  (D) . 


2 . 3  Relaxation  (PRn) 

Theorems  1  and  3  are  useful  in  providing  a  choice  of  Lagrange 

multipliers  as  a  starting  point  in  solving  a  relaxation  of  problem  (P) 

at  the  root  node.  Theorem  2,  similar  to  Theorem  1,  provides  results  for 

a  partial  solution  of  problem  (P),  i.e.,  at  a  node  other  than  the  root 

node  where  some  of  th.e  x..  and  y,  have  been  fixed  at  1  or  0  . 

ij 

Theorem  1  is  important  in  pointing  out  that  .i  certain  set  of  Lagrange 
multipliers  v  such  tliat  v.  >  b,  for  all  k  would  provide  an  optimal 

choice.  Theorem  3  narrows  this  choice  to  v  =  b  for  all  k  for  a 

k  c 

specific  situation,  i.e,,  when  the  resulting  solution  is  feasible  for 
problem  (P) . 

Letting  Vj^  =  bj^  for  all  k  ,  problem  (LR^)  becomes: 

f 


Minimize 

Z  I  c  .  . 

i  j 

X  .  . 

ij 

(28) 

subject  to 

V3 

(2) 

X.  .  =  0 
ij 

or  1 

Vi,j 

(6a) 

where  c  .  . 

=  a  .  .  +  Z 

12  V 

b.  r . . 

K  ijk 

(29) 

Note  that  problem  (LR^^)  is  very 
the  sum  of  the  minimum  (over  i) 


easy  to  solve;  its  optimal  value  is  just 
for  all  j  ,  i.e. , 


Z(LR  )  =  Z  min  {c . 


(30) 


We  solve  this  problem  as  a  starting  point  at  the  root  node  in  our  branch- 
and-bound  procedure. 
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so  is  fixed  at  1  and  kcT  if  =  1 


To  return  to  the  relaxation,  problem  (LR£  ,  we  must  make  a  choice 

of  the  vector  v  of  Lagrange  multipliers.  Of  course,  we  would  like  to 
use  an  optimal  choice,  i.e.,  a  vector  v  that  solves  problem  (3^)  • 

Recall,  however,  that  Theorem  2  did  not  provide  us  any  useful  information 
about  tlie  optimal  value  of  v^^  if  keT  and  y^^  “  1  .  To  simplify  our 

approacn  and  have  recourse  to  the  results  of  Theorems  1  and  3,  we  choose 

V  =  0  if  keT  and  y  =  1  .  Note  that  there  a  re  no  kcT  such  that 

k  k 

V  =  0  because  of  practical  considerations  and  because  our  branching  rule 
-'k 

only  results  in  fixing  values  at  1  .  Problem  (LR^  now  takes  the 


Minimi  ;e  IIx..  (  a,.  +  Z  v,  r..,\-Z  y,  /v,  -  b,  I 

i  j  kef  kei  4  ^  7  kGT 


h^(35> 


(..R£^-)  i  subject  CO  Z  x.j  =  1  , 


Xfj,  yj^  =  0  or  1  for  all  (i,j)eS  ,  keT  . 


Note  that  in  this  problem  (LR.  -) ,  v  ■=  0  if  kcT  .  .\lso  note  how  closely 

X/  f  V  K 

it  resembles  problem  (LR^),  the  relaxation  at  the  root  node.  As  in  that 

case,  wo  would  like  the  lower  bound  Z(LR.  -)  to  be  as  large  as  possible, 

X ,  v 


i.e.,  we  seek  v  to 


Maximize  ^Z(LR£  -)J 


V  >  0 


Because  of  the  close  similarity  of  problems  (LR^  -)  and  (LR^), 

it  is  possible  to  obtain  results  about  problem  that  are  analogous 

to  those  obtained  about  problem  (D) .  We  state  these  results  as  Theorems 
4  and  5.  Their  proofs  are  omitted  because  they  follow  precisely  the 


lu'ool  s  ol  I’l'.uoroi'.is  1  anil  1,  respectively,  and  their  v.ilidity  follows 
from  the  fact  that  problem  (i,R  -)  is  essentially  the  same  as  problem 

,  V 

(i.Ky)  but  involves  only  the  free  variables. 

H'-i’l-'Jll'ILA-  i'here  exists  an  optimal  solution  to  problem 

in  which  v  '•  b,  Lor  all  keT  . 

K.  --  k 

'k 

Theorem  5:  Let  (X  ,  Y  )  solve  problem  (LR„  -)  for 

A/  ,  V 

_  k  ie 

for  all  keT  .  If  (X  ,  Y  )  satisfies 

(5')  for  all  kcT  ,  there  exists  an  optimal  solution 

to  problem  (Dt)  in  which  v,  =  b  for  all  keT  . 

■It  k  k 

Just  as  Theorems  1  and  3  motivated  us  to  use  the  relaxation 
problem  obtain  our  lower  bound  at  node  1  ,  Theorems  A  and  5 

motivate  us  to  set  foi-  relaxation  problem  (LR^  -) 

to  obtain  our  lower  bound  at  node  Z  .  With  this  specification,  problem 
(LR  -)  becomes 

X, ,  V 


Minimize  .1'  L  c...  x  +  FC„ 

.  .  ijL  ij  i 

1  J 


(PR,)  -I  subject  to  L  X.,  =  1 


X .  .  =  0  or  1 
ij 


V(i,j)eS  , 


(23a) 


where 


c  ,  .  0  =  a  .  .  +  L  b,  r  .  . , 

iji(  ij  k  ijk 


*'  k-i  \ 


1  -  6 


kij  i  i  1; 


and  the  fixed  cost  is  given  by 


FC„  =  I  b ,  =  I  6,  0  b.  . 
^  keT  k=l 


(39) 
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This  specific  relaxation,  problem  (PR^) ,  is  oi  the  same  form  as 
problem  (LR,  )  and  is  equally  easy  to  solve  in  one  pass.  Its  optimal 
value  Z(PRj^)  serves  as  the  lower  bound  at  node  P.  .  Note  that  for  il=l 
problem  fhe  same  as  problem  (LR^^)  . 

It  is  clear  that  setting  each  Lagrange  multiplier  v^^  to  b^^ 

for  keT  and  to  0  for  keT  is  not  generally  optimal  in  terms  of 
achieving  the  tightest  lower  bound  (except  as  per  Theorem  3) .  But  it 
provides  a  good  starting  point  in  seeking  an  optimal  vector  v  and  it 
provides  an  easily  calculated  lower  bound  at  each  node  of  our  branch- 
and-bound  procedure.  The  question  of  hov/  to  improve  upon  this  choice 
of  multiplier  values  will  be  discussed  in  Chapter  6. 


3 .  MET’.'.ODOLOGY  FR.\MEWORK 


The  branch-and-bound  procedure/metliodo] ogy  developed  to  solve 
problem  (P)  uses  Lagrangian  relaxation  (f’R£)  a  basic  step.  The 

branching  rule  dictates  which  x,.  variable  to  branch  on  at  each  node, 

IJ 

In  addition,  there  are  certain  rules  (e.g.,  the  capacity  rule  and  the 
bounding  rule)  which  contribute,  significantly,  in  improving  the  overall 
efficiency  of  the  procedure. 


Some  basic  terras  such  as  fixed  and  free  vaiiablcs,  partial  solution 
and  its  completion  were  introduced  in  the  previous  chapter.  This  chapter 
first  provides  a  preliminary  discussion  of  the  branch-and-bound  methodology, 
[Ceol’frion  (1967),  and  Geoff rion  and  Ilarsten  (1972)].  Representation  and 
storage  of  the  x_  variables  for  branching  and  backtracking  is  described 

in  order  to  pro\ ide  continuity  and  consistency  with  the  computer  program 
covered  in  Chapter  i.  'This  is  followed  by  a  description  of  the  major 
components  of  the  branch-and-bound  methodology. 


Branching  and  backtracking  is 

branching  commences  by  tixinp,  tlie  x 

rule)  to  1  and  moying  to  thi-  left 
1  ix  the  corresponding  x..  vari.ible 


done  on  the  x.,  variables.  The 
ij 

variable  (selected  by  the  brandling 

branch  node.  When  backtracking,  we 
at  0  and  move  to  the  right  branch 


node  (if  tlie  rignt  branch  node  has  not  already  been  explored).  An 

variable  can  also  be  fixed  at  0  or  1  by  rules  other  than  the  branching 
rule.  The  capacity  rule  and  tlie  bounding  rule  arc  two  such  rules  employed 
in  our  methodology. 


Figure  4a  shows  a  branch-and-bound  tree.  The  variables  fixed 

at  0  or  1  at  any  node  due  to  rules  other  than  the  brancl.ing  rule  are 
shown  in  parenthesis  at  the  appropriate  node. 
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yode  ©  is  the  root  node  and  also  the  par.mt  node  for  nodes 
0  and  © 

ilode  @  is  the  parent  node  for  nodes  ©  and  ©.  etc. 


Figure  4a. 

A  branch-and -bound  tree  illustration 


Node  (f/) 

1 

2 

3 

4 

5 

6 
7 


Partial  .Solution  (Sj^) 


<P 

{103,  -  301,  -  ^01} 

{103,  -  301,  -  ■^Ol,  402} 

{103,  -  301,  -  401,  -  A02,, 

{103,  -  301,  -  ^01,  -  402, 

{103,  -  301,  -  401,  -  402, 

{-  103} 


201,  101} 

201,  101,  204} 
201,  101,-204} 


Partial  solutions  for 


Figure  4b. 

the  above  illustration  (Figure  4a) 
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For  problem  (PR^)  .  a  partial  solution  corresponding  to  set  S  at 

node  I  ,  i.e  ,  S,,  contains  x..  variables  assigned  values  of  1  or 

i.  ij 

0  .  For  simplicity  in  the  computer  program,  an  variable  fixed  at 

1  is  represented  as  (100  i  +  j)  ,  whereas  an  variable  fixed  at 

0  as  -  (100  1  +  j)  ,  e.g.,  x^,,  =  1  and  x^2  "  ^  represented  as 

302  and  -  302  respectively.  Since  branching  is  done  on  x^^  variables, 

it  is  necessary  to  make  a  distinction  between  x^^  variables  fixed  at  1 

due  to  the  branching  rule  and  those  fixed  at  1  due  to  the  other  rules. 

We  make  this  distinction  by  underlining  the  positive  number  to  represent 

an  X  .  fixed  at  1  due  to  the  ocher  rules.  For  example,  204,  -  301, 
ij 

1^3  represent,  respectively,  X2^  =  1  due  to  the  branching  rule,  x^^  "  0 
due  to  the  branching  rule  or  any  other  rule,  and  x^^  ~  1  due  to  a  rule 
other  than  the  branching  rule. 


Figure  4h  sl\ows  the  partial  solutions  of  tb.e  branch-and-bound 

tree  in  Figure  4a. 

Implicit  enumeratLon  involves  generating  a  seciuence  of  partial 
solutions  and  simultaneously  considering  all  completions  of  each.  For 
our  minimization  problem,  we  start  with  an  initial  solution  having  a 
very  large  value  (infinity)  as  an  initial  upper  bound.  As  the  computations 
proceed,  feasible  solutions  (those  satisfying  the  capacity  constraints) 
are  disccjvered  from  time  to  time,  and  tlie  best  one  yet  found  is  retained 
as  an  incumbent  solution  with  the  corresponding  value  as  the  best  upper 
^omid.  It  may  happen  that  for  a  given  partial  solution  we  can 

determine  a  best  completion  of  .S^  ,  i.e.,  a  feasible  completion  that 

minimizes  the  objective  function  value  among  all  feasible  completions 
of  .  If  such  a  best  feasible  completion  is  better  than  the  best 

upper  bound,  then  it  replaces  the  l.ittei .  Or  wo  may  be  able  to  determine 
that  has  no  feasible  completion  bettor  than  the  incumbent.  In  either 

case,  we  can  f .1  tiioin  .  (V'arious  situations  of  1  athoming  and  back- 
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tracking  in  our  branch-and-bound  procedure  are  doacribed  in  the 
following  discussion.)  All  completions  of  a  fatho-ned  partial  solution 
have  been  implicitly  enumerated  in  the  sense  that  they  can  be 

excluded  from  further  consideration  (with  the  exception  of  the  relevant 
best  feasible  solution  of  if  it  has  been  retained  as  the  best 

upper  bound) . 

In  our  branch-and-bound  procedure,  at  any  given  node  where  we  can 
fathom  ,  we  backtrack  to  the  parent  node  and  move  to  the  right-hand 
branch  (if  that  branch  has  not  already  been  explored)  by  fixing  the 
appropriate  variable  at  0  .  However,  if  the  right-hand  branch 

has  already  been  explored,  we  continue  backtracking  to  a  parent  node  where 
we  can  move  to  a  right-hand  branch.  For  example,  in  Figure  4a,  when 
backtracking  from  node  3,  we  move  to  the  parent  rode  2,  and  to  the  right 
to  node  4  by  setting  x^2  “  ^  .  However,  when  backtracking  from  node  6, 

we  move  back  to  node  4,  then  back  to  node  2,  then  back  to  node  1,  and  to 
the  right  to  node  7  by  setting  =  0  . 


On  the  other  hand,  if  the  partial  solution  £,  cannot  be  fathomed, 

X, 

we  branch  to  the  left  and  augment  S.  by  fixing  a  free  variable  x  at 

1  (based  on  the  branching  rule),  and  then  we  try  to  fathom  the  resulting 

partial  solution.  In  addition  to  the  one  variable,  selected  by  the 

branching  rule,  some  other  free  x. ,  variables  can  also  be  fixed  at 

ij 

0  or  1  according  to  the  application  of  rules  other  than  the  branching 
rule.  Note  that  this  can  also  happen  when  backtracking^,  i.e.,  when 


ha.s  been  fathomed  and  we  backtrack  and  move  to  the  right  by  setting 

the  appropriate  x..  variable  to  0  . 
tj 


Let  us  consider  examples  of  both  situations,  i.e.,  when 
has  not  been  fathomed  and  when  has  been  fathomed.  In  Figure  4a 

we  cannot  fathom  (i.e.,  S  at  node  1),  so  we  move  to  node  2  by 
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.mii'.uenL in;;  S,  by  I  ixinj^  =  1  based  on  the  branching  rule,  and 

by  fixing  x^j^  =  0  and  x^^  =  0  based  on  the  aapjication  of  tlie  other 

rules.  Similarly,  we  move  from  node  2  to  node  3  by  augmenting  S2  by 

tixiag  x^,  =  I  .  As  an  example  of  backtracking,  wher  we  fathom  , 

we  move  back  to  tlie  parent  node  2,  and  to  the  right  tcj  node  4,  getting 
a  new  partial  solution  S,  by  re[)lacing  x,  =  1  with  x ,  ~  =  0  ,  and 

further  augm.Miting  it  by  fixing  x,,,  =  0  and  x,,  -  1  based  on  the 

iX 

application  of  the  otiier  rules. 

if'mputation.iliy ,  the  storage  and  update  of  partial  solution 
is  easily  accompJished  by  considering  Figure  4b.  If,  at  a  given  node, 
the  partial  solution  has  not  been  fathomed,  e.g.,  at  node  4, 

determine  the  next  branching  variable  by  using  the  branching  rule, 
i.e.,  ,  and  augment  by  adding  204  as  the  last  entry.  Also, 

augment  S^^  w- th  any  other  free  x^^  variables,  if  appropriate, 

depending  on  the  application  of  the  other  rules.  Now,  consider  the 
case  where  the  partial  solution  has  been  fathomed,  e.g.,  at  node  6, 

and  we  backtrack;  starting  with  the  last  entry  in  S^,  .  we  consider  one 

entry  at  a  time,  going  backwards,  until  wc  find  a  posi  ive  number  which 
is  not  untlcrlinod.  In  our  example,  it  is  103.  In  otlii  r  words,  we  must 
branch  to  ti‘2  right  by  fixing  =  0  ,  i.e.,  vie  replace  103  with  -  103 

and  we  are  at  node  7.  Should  we  find  that  we  have  no  positive  number, 
the  procedure  terminates  since  we  are  b.ick  at  the  root  node  and  the 
rig.tii  branch  has  already  been  explored.  This  happens  when  backtracking 
from  node  7. 

In  t  lie  br.inc h-and-bound  procedure  wc  generate  a  sequence  of 
partial  solutions  as  we  move  from  one  node  to  another.  This  sequence  is 
iioii-i  cduud.int  in  the  sense  that  no  completion  of  a  partial  solution 
ovor  duplicates  a  completion  of  a  previous  parti.il  solution  that  has 


been  fa  tliomed. 
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Since  one  of  the  x  values,  for  each  j  ,  musf  be  1  ,  a  total  of 

ij 

(2m-l)  nodes  are  theoretically  possible  for  complete  enumeration. 
However,  most  of  the  solutions  may  be  infeasiblt  because  of  the  capacity 
constraints.  The  branch-and-bound  procedure,  thro-.gh  a  judicious  choice 
of  branching  variables,  and  elimination  of  certain  infeasible  and  non- 
optimul  assignments  tlirnugh  various  rules,  turns  out  to  be  a  practical 
and  computationally  efficent  algorithm.  The  various  components  of  this 
procedure  are  described  next.  Detailed  procedural  steps  and  the  solution 
of  a  test  problem  will  be  covered  in  Chapter  4. 

3 . 1  Bounds 

3.1.1  Lower  Bound 

At  a  given  node  Z  in  the  branch-and-bound  tree,  a  lower  bound 
(LOWB)  is  obtained  by  solving  relaxed  problem  (PR^)  • 


LOWB  =  Z(PRJ 

£ 


(40) 


Recall  that  problem  (PR^^)  is  very  easy  to  solve  by  considering  the 
minimum  c...  over  chose  j's  for  which  x,.  is  n^t  fixed  at  1  , 
i.e.,  jeW  ,  where  W  is  the  complement  of  W  defined  by  expression  (32), 


Z(PR  )  =  I  c..„  +  E 
*  JCU  JEW 


(i,  j)es 


(41) 


where  given  by  expression  (38),  I.e., 


"ijSL  °  ^j  \  ’^ijk  ’ 


(38) 


and  the  fixed  cost  (PC^)  is  given  by  expression  (39),  i.e.. 


"S  -  “la  \  • 

k 


(39) 


where  is  given  by  expression  (34). 


i 


I 


Note  that  it  none  of  the  x.  variables  is  fixed  at  1  ,  as  is 

ij 

i’onerally  the  case  at  the  root  node,  then  all  6,  ,  =  0  ,  and,  therefore, 

kj. 

P 

FC  -  0  ,  and  c .  -  a . .  +  Z  b  r .  .  Z (PR  )  is ,  then ,  simpl v  the 

ijl  ij  k  ijk  1  .  t  - 

middle  part  of  expression  (Al)  .  Vic  use  the  term  "generally''  because 
it  is  possible  that  tlie  capacity  rule  could  force  certain  x 

ij 

variables  to  1  (or  0)  at  the  root  node,  prior  to  solv'ng  the  relaxed 
problem  (PR^)  . 

1.1.2  Upper  Bou.nd 

At  any  given  node  i  ,  let  X  =  represent  the  solution  of 

problem  (FR^)  .  ft  this  solution  is  feasible  for  problem  (P),  i.e., 
if  X  satisfies  the  capacity  constraints  (5)  or  (5') 

I  \  \  Vk  ,  (42) 

Xi.cX 

where  y,  =  i  if  Z  d...  x  >  0  , 

k  .  j  ijk  ij 

x..t:X  (4  3) 

=  0  otlu-rwise  , 


tlien  the  value  of  problem  (P)  corresponding  to  this  solution 
gives  an  upper  bound  (V'PB)  : 


r  L  a.  .  X.  .  +  I  b  V,  , 
•  •  ij  ij  ,  k  k 

1  J  k 

X .  .ex 


wliere  y^^  is  defined  by  (AJ)  . 


3.1.1  i’est  Upper  Ihnird 


A  current  Lo’..i'st  upper  bound  is  retained  .as  tlie  best  upper  bound 
(bllB),  the  corresponding  solution  X  representing  the  incumbent  solution. 
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The  brancl'.-and-bound  procedure  is  initiated  by  assumz'ng  a  very  large 
value  as  Che  best  upper  bound,  and  is  replaced  by  better  (lower)  values 
as  the  procedure  continues. 

A  positive  fractional  value  c  can  be  specified  if  a  sub-optimal 
solution  is  acceptable.  For  example,  for  e  =  0.001  ,  the  resulting 
solution  value  is  guaranteed  to  be  within  0.1  percent  of  the  optimal 
solution  value.  When  e  is  non-zero,  the  adjusted  best  upper  bound  (BUBS) 
is  defined  as: 

BUBS  =  BUB/(1  +  e)  .  (45) 

Obviously  when  e  =  0  ,  BUBS  =  BUB  . 


3.2  Facility  Usage  Rule 

This  rule  is  used  to  identify  facilities  forced  into  usage  at  a 
given  node  2.  and  hence  fix  corresponding  free  variables  y^^  at  1  . 


For  a  partial  solution  ,  define 


^jk2  =  ‘^ijk  ’ 


=  min  d . . ,  if  j  eW  . 

(i,j)eS 


(46) 


The  facility  usage  rule  states 
is  not  already  fixed  at  1  ,  if 


that 
I  d 
j 


jk2 


for  any  facility  k  , 
>  0  ,  tiien  facility 


where 
k  is 


forced  into  usage  and,  therefore. 


should  be  fixed  at 


1  . 


This  rule  is  applied  at  every  node  prior  to  applying  the  capacity 
rule.  In  other  words,  this  rule  is  applicable  to  capacitated  as  well  as 
uncapacitated  problems. 

3 . 3  Capacity  Rule 

This  rule  is  designed  to  "exclude"  infeasible  assignments  prior 
to  solving  the  relaxed  problem  (PR£)  •  This  is  done  by  exploiting  the 


-  J 


rei.iL  iuiKshi]'  beLweeii  the  capacities  requited  (d..,)  and  the  capacities 

1  j  k 

available  )  ier  a  ^iven  partial  solution  ol  problem  (P). 

The  capacity  rule  states  that  for  a  facility  k.  and  an  activity 
j  ."exclude" a  free  x..  variable  (i.e.,  fix  it  at  0)  for  which 

ij 


(i, j)f  S 


f47) 


where  d.,  .  is  defined  by  expression  (46).  The  right-hand  side  of  this 
J 

inequality  (47),  when  positive,  represents  the  available  capacity  at 
facility  k  .  Tiie  left-hand  side  shows,  for  a  given  j  ,  the  difference 
between  a  d...  corresponding  to  a  free  x.,  variable  and  d,,„  . 

ijk  c  to 

If,  fur  a  specific  d...  ,  this  difference  is  more  than  the  available 

ijk 

capacity,  the  corresponding  free  x,.  variable,  if  fixed  at  1  ,  would 

result  in  an  infeasible  solution.  Thus,  by  looking  ahead,  we  can  exclude 
sucii  a  free  variable  by  assigning  it  a  value  of  0  . 

Note  that  if  the  right-lrand  side  of  expression  (47)  is  negative, 
then  any  completion  of  sucli  a  partial  solution  will  be  infeasible  and 
we  backtrack  in  our  branch-and-bound  procedure. 

The  capacity  rule  is  applied  to  all  the  facilities  by  considering 
one  facility  at  a  time.  The  cycle  of  examining  all  tht  facilities 
continues  until  no  more  assignments  can  be  excluded.  During  the  course 
of  application  of  this  rule,  if  all  but  one  of  the  free  x^^  variables 

h.ive  been  excluded  (fixed  at  0)  for  a  given  j  ,  then  that  particular 

X..  variable  is  fixed  at  1  because  of  constraints  (2),  i.e.,  each 
J-l 

activity  j  must  be  assigned  to  one  and  only  one  design  1  .  The 
partial  solution  is  updated  accordingly  to  reflect  the  x^^  variables 

fixed  at  0  or  i  due  to  the  application  of  the  capacity  rule. 


The  capacity  constraints  for  an  uncapacitated  problem  are  not 

Hence,  the  capacity  rule  is  useful  only  for  capacitated  problems. 


active. 
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3 . 4  Branch ing  Rule 

This  rule  provides  the  choice  of  the  variables  on  which  to 

branch.  If  the  partial  solution  at  a  given  node  H  is  not  fathomed,  we 
branch  further  by  fixing  a  free  variable  at  1  and  moving  to  the 

left  branch  node. 

According  to  the  branching  rule  the  choice  of  the  branching 

variable  depends  on  the  c  .  values  and  is  such  that  the  corresponding 

ij  X. 

x^j  ,  if  perturbed,  has  the  maximum  impact  on  the  optimal  value  of 


problem  (PR£)  • 


For  a  given  j  ,  define  c  ,  the  minimum  permissible  c .  .  j.  , 

ijjx.  ijx, 


the 

second 

smallest 

permissible 

‘"iU  ’ 

=  min 
i 

‘'ij£ 

for 

jeW  and 

(i.:)es 

(48) 

=  min 
i 

for 

j  eW  and 

(i,j)es 

(49) 

For  each  jeW  ,  define  D.„=c.  ..  -c.  ..  .  (50) 

jx.  1-2^^  ^2.^^ 

Our  branching  rule  states  that  a  free  x^^  variable  corresponding  to 
c.  .g  such  that  D.-  is  maximized  over  all  j  ,  is  selected  as  the 

J  X, 

next  branching  variable  and  assigned  a  value  of  1. 

3. 5  Bounding  Rule 

This  rule  is  designed  to  "exclude"  certain  non-optlmal  assignments. 
These  assignments  cannot  lead  to  an  optimal  solution  as  we  branch  from 
one  node  to  the  next  left  branch  node. 


The  bounding  rule  states  that  a  free  variable  should  be 

excluded  (by  assigning  it  the  value  0)  for  which 


(c  .  .  .  -  c  .  .„)  >  (B’JHS  -  LOWB)  for  jcW  and  (i,j)cS 


(51) 


ij'- 

L 

where  c  .  ,  BL'BS,  and  LOWB  are  given  by  exi'.ressions  (48),  (45), 
and  (40),  respectively. 


Thus,  by  lotiking  ahead,  we  exclude  those  .assignments  which  will 
provide  lower  bounds  higher  tlian  BUBS. 

The  bounding  rule  is  applied  to  each  jeW  just  prior  to  selecting 
tlio  X  varial)ie  for  branching  to  the  left. 

Ij 


As  in 
in  e.xcluding 

a  given  jCW  , 


the  case  of  the  capacity  rule,  if  Che  bounding  rule  results 

(fixing  at  0)  all  but  one  of  Che  free  x..  variables  for 

ij 

then  that  particular  x  .  variable  is  fixed  at  1  .  Also 

1.1 


the  partial  solution  is  updated  accordingly  to  reflect  the  x.. 


variables 


fixed  at  0  or  1  due  to  the  application  of  the  bounding  rule. 


^ ^  B^t^ kt racking  Rules 


If  a  partial  solution  at  a  given  node  has  been  fathomed,  we  back¬ 
track.  The  backtracking  rules  are  typical  of  a  branch-and-bound  procedure. 
In  addition,  the  application  of  the  capacity  rule  and  the  bounding  rule 
can  lead  to  backtracking.  The  criteria  for  backtracking  include  the 
fo  1.  lowing. 


(a)  When  applying  the  capacity  rule,  if  the  available 

cajiacity  given  by  the  right-hand  side  of  inequality 
(47)  is  negative,  i.e.,  (.s^^  -  I  5^^^^)  <  0  , 

then  backtrack.  ^ 


(b)  If  LOWB  ■  HUBS,  then  backtrack.  Otherwise  compute 
Ui'B  if  the  solution  is  feasible  in  problem  (R). 

Then  update  BUR  and  BUBS  if  Ul’B  <  BUB;  and  back¬ 
track  if  LOWB  =  BUBS. 

(c)  If  furtlicr  brandling  Is  not  possible,  then  backtrack. 

This  c.in  happen  due  to  the  cap.icity  rule,  the  bounding 

rule,  or  the  hr.inching  rule  if  the  updated  partial 

solution  is  such  th.it  no  further  brandling  is  possible, 

i.e.,  X,,  variables  are  fixed  at  1  for  all  j  ,  or 

tj 

equ  Iva  lent  Ly  ,  Vi  =  <ti 


Wl'.ea  ony  ot  the  backtracking  criteria  apply,  we  backtrack  to 
the  parent  node  and  move  to  the  right  branch  node  (if  the  right  branch 
has  not  already  been  explored)  by  fixing  the  appropriate  x^^  variable 

at  0  .  If  the  right  branch  has  already  been  explored,  we  continue  back¬ 
tracking  to  a  parent  node  where  we  can  move  to  a  right  branch  node. 

The  branch-and-bound  procedure  terminates  when  we  backtrack  to  the 
root  node  and  find  that  the  right  branch  node  has  already  been  explored. 


f 


r 


A.  COMPUTATIONAL  STEPS  AND  THE  COMPUTER  PROGRAM 


A  computer  program  called  ZIPCAF  (an  acron^TH  for  Zero-one  Integer 
Program  for  raultiac Liv ity  multifacility  Capac ity-i  onstrnined  Assignment 
Problems)  implementing  tiie  brancli-and-bound  methodology  has  been  developed. 

Detailed  procedural  steps  and  guidelines  to  use  the  computer 
program  are  described  in  a  separate  document  [Chhabra  and  Soland  (1980)] 
titled  "Program  Description  and  User's  Guide  for  ZIPCAP — a  Zero-one 
Integer  Program  to  solve  multiactivity  muitifncillty  Capac Ity-constralned 
Assignment  Problems."  Specifically,  the  document  includes: 

.  Problem  fornul.ition  (]’)  and  potential  areas  of 
appl icat ion 

.  Overall  flow  diagram  and  detailed  procedural  steps 
for  the  computer  program 

.  Program  listing  and  dictionary  of  t'.ie  symbolic  names. 

The  listing  includes  extensive  use  of  comment  cards 
to  explain  various  computational  steps. 

.  User  information  including 

schematic  diagram  of  the  deck  structure, 

detailed  instructions  for  the  job  control  (JCL) 
cards,  program  parameter  card,  program  options 
card,  and  the  various  other  input  data  cards. 

.  Three  test  problems  to  demonstrate  the  use  of  the 
program.  The  display  Includes  coded  input  and 
annotated  outputs  reflecting  the  use  of  selected 
prop7  ram  o p t  ions  . 

As  mentioned  ear  ler,  ZIPCAP  is  primarily  designed  for  capacitated 
problems.  However,  uncapacitated  problems  can  be  solved  as  a  special 
case,  and  this  is  demonstrated  by  including  an  uncapacitated  test 
problem. 


A7 
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Because  of  the  extensive  coverage  of  the  program  description  and 
user  guidelines  in  the  above  document,  this  chapter  provides  only  an  over¬ 
view  of  the  computer  program,  including  an  overall  flow  diagram,  and  a 
summary  of  the  program  options,  in  order  to  provide  continuity  in  this 
document.  In  addition,  a  step-by-step  description  oi  a  test  problem  is 
presented  to  demonstrate  the  use  of  the  various  components  of  the  branch- 
and-bound  methodology.  The  computer  printout  showing  step-by-step  details 
is  obtained  by  use  of  one  of  the  program  options.  The  use  of  this  option 
to  display  detailed  steps  in  this  document,  in  fact,  complements  the  use 
of  the  various  options  demonstrated  in  the  other  document. 

4 . 1  The  Program 

Figure  5  presents  a  simplified  flow  diagram  of  the  branch-and- 
bound  procedure.  The  major  computational  steps  for  the  computer  program 
are  numbered  in  circles.  These  steps  are  essentially  based  on  the 
methodology  components  described  in  the  previous  chapter.  A  step-by-step 
description  has  been  Included  in  the  other  document  [Chhabra  and  Soland 
(1980)1. 

The  computer  program  ZIPCAP  is  written  in  FORTRAN  IV,  and  has 
been  developed  and  tested  on  an  IBM  3031  at  the  George  Washington 
University.  The  program,  comprising  about  480  lines  is  currently 
dimensioned  for  a  maximum  problem  size  of  35  designs  (m/ ,  35  activities 
(n)  and  30  facilities  (p) .  The  program  size  to  execute  a  problem  has 
two  components:  one,  due  to  the  program  Itself,  comprising  173  K  bytes, 
and  the  other  aependent  on  the  dimensions  of  the  arrays  given  by  the 
following  functional  relationship. 

f(m,.i,p)  =  4[(p+4)mn  +  (m+5)p+9n]  bytes 

The  computer  program  listed  in  the  other  document  has  since  been 
further  improved.  The  basic  improvement  has  been  the  addition  of  the 
facility  usage  rule.  This  rule,  as  described  in  Chapter  3,  is  applied 
both  to  capacitated  and  uncapacitated  problems  just  before  the  applica¬ 
tion  of  the  capacity  rule.  For  completeness  of  this  document,  a  revised 
program  listing  is  included  in  Appendix  A.  It  may  le  mentioned  that 
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Figure  5. 

Simplified  flow  diagram  for  the  branch-and -bound  procedure 
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the  revised  program  solves  the  test  problems  incl'ided  in  the  other  document 
more  efficiently  —  in  less  time  and  in  fewer  noc  es  (with  an  average 
reduction  in  nodes  of  31  percent).  The  improvement  in  efficiency  seems 
to  result  from  the  "multiplicative"  effect  of  the  various  rules.  Another 
improvement  made  is  that  the  computer  printout  always  displays  th*j  node 
number  (IBNOD)  at  which  the  best  upper  bound  changes  (improves)  and  the 
corresponding  values  of  the  best  upper  bound  (BUB)  and  the  adjusted  best 
upper  bound  (BLBS). 

ZIPCAP  provides  numerous  options  to  the  program  user.  These  options, 
described  in  the  other  document,  are  summarized  in  Table  3. 

Option  ICAPR,  the  capacity  rule,  is  automatically  skipped  by  the 
program  when  solving  an  uncapacitated  problem.  Option  ISTEP,  the 
intermediate  steps'  listing,  even  when  skipped,  provides  information  on 
the  total  number  of  nodes  explored.  A  summary  listing  provides  necessary 
information  tc  construct  the  branch-and-bound  tree,  whereas  a  detailed 
listing  of  the  intermediate  steps  is  useful  when  changing  or  debugging 
the  program. 

Option  EPS,  the  optlmal/suboptimal  solution,  provides  the 
flexibility  of  obtaining  a  suboptlmal  value  guaranteed  to  be  within 
a  specified  fraction  of  the  optimal  value.  The  resulting  solution  may 
he  suboptlmal  but  could  provide  a  considerable  saving  in  terms  of 
exploring  fewer  nodes  in  comparison  to  those  necessary  for  obtaining 
an  optimal  sclutj.on. 

Option  ET,  by  providing  important  information  at  a  specified 
elapsed  time,  is  useful  in  a  situation  where  the  t>ital  time  allocated 
to  solve  a  problem  may  not  be  sufficient  and  the  program  terminates 
before  verifying  an  optimal  solution.  The  information  provided  by 
this  option  includes  an  updated  partial  solution  showing  the  x^^ 

variables  fixed  at  0  or  1,  at  the  current  node  being  explored  at 
the  specified  tire  ET.  By  looking  at  the  first  few  variables  displayed 
in  the  partial  solution  of  the  current  node,  it  is  possible  to  assess 
the  extent  of  the  branch-and-bound  tree  explored  until  time  ET.  For 


iriLTa  >1  Ti  ni"^r 


SUMMARY  OF  ZIPCAP  OPTIONS 


Skip  this  option 


example,  in  view  of  the  terminology  in  Figure  4b  (Chapter  3),  if,  at  an 

arbitrary  node,  the  first  term  of  the  partial  solution  is  positive, 

i.e.,  the  x  .  variable  has  value  1,  then  we  are  still  in  the  left 
ij 

half  of  the  total  branch-and-bound  tree.  If  the  first  term  is  negative, 
i.e.,  the  x^j  variable  has  value  0,  then  we  ar'^.  in  the  right  half  of 

the  total  branch-and-bound  tree  and  have  explored  half  of  the  total 
(theoretical)  solutions  corresponding  to  the  left  half  of  the  tree. 

If  the  first  two  terms  are  negative,  i.e.,  the  fi'-st  two  x^^ 

variables  have  value  0,  then  one  quarter  of  the  total  (theoretical) 
solutions  remain  to  be  explored,  since  we  are  in  the  next  right  half  of 
the  right  half  of  the  total  branch-and-bound  tree,  as  illustrated  in 
Figure  6. 


Each  j  egment  represents  1/4 
of  the  total  (theoretical) 
solutions 

Each  segment  represents  1/2 
Left  half  of  Right  half  of  of  total  (theoretical) 

the  "total"  the  "total"  solutions 

branch-and-  branch-and- 

bound  tree  bound  tree 


Figure  6 

Illustration  for  estimating  the  extent  of  the  branch' 
and-bound  tree  explored 
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Ki-'call  froH’  Chapter  3,  that  a  total  of  (2m-l)''  nodot  are  theoretically 
possible.  Tlurs,  it  the  first  g  [g  <  (m-l)n]  terois  at  an  arbitrary 

node  are  negative,  then  theoretically  about  [ (2n— l)*'/^g]  nodes  remain 
to  be  explored. 


4 . 2  An  Illustrative  Example 

We  consider  a  capacitated  test  problem  with  five  designs  (m) , 
four  activities  (n) ,  and  eight  facilities  (p)  to  demonstrate  the  use 
of  the  branch-and-bound  procedure  and  the  computer  program. 

The  computer  printout  for  this  problem  shoving  step-by-step 
details  for  a  couple  of  nodes  is  presented  in  Appendix  B. 

As  shown  in  the  beginning  of  the  printout,  the  options  selected 

are: 


.  lINPT  =  1,  i.e.,  list  the  input  data 

.  ICAPR  =  1,  i.e.,  use  Che  capacity  rule 

.  ISTEP  =  2,  i.e.,  list  detailed  intermediate  steps 

.  lUNCAP  =  1,  since  this  is  a  capacitated  problem 

.  EPS  =0.0  implying  that  an  optimal  solution  is  desired 

.  ET  =  0.0  since  a  detailed  listing  of  intermediate  steps 

will  be  available. 


Kol lowing  the  listing  of  the  options,  input  data  listed  for  the  problem 
include  variable  costs  a..  ,  fixed  costs  b,  ,  available  capacities 

,  and  cap.'cities  required  values  ire  generated 

by  thi-  computer  program. 


The  computer  program  follows  the  procedural  steps  marked  in  the 
flow  diagram  presented  in  Figure  5.  These  steps,  ilong  with  the 
relevant  terminology  used  in  the  computer  printout,  are  described  below 
for  a  couple  of  node.s,  followed  by  a  complete  br.inch-ana-bound  tree  for 
this  problem.  As  mentioned  earlier,  a  dictionary  of  the  symbolic  names 
usi'd  in  the  computer  program  is  included  in  the  other  document. 
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Node  1 

Seep  J:  Initialize. 

Initialize  BUB  =  9999999.0,  and  since  EPS  =  0.0,  BUBS  =  BUB. 

Also  S  =  0  and  W  =  .  In  the  computer  printout,  vector 

FIX(J)  represents  the  sCi.  W,  and  matrix  CX(I,J)  ;epresents  both,  fixed 

and  free  x  .  variables.  In  the  CX(I,J)  matrix,  an  x.,  variable  fixed 
-LJ  iJ 

at  1  or  0  is  represented  as  1  or  2,  respectively,  and  a  free  x^^  variable 

is  represented  by  the  value  0.  Initially,  all  the  variables  are  free 

as  sho^^7n  by  n.atrix  CX(I,J)  in  the  printout. 

Step  2:  Apply  the  facility  usage  rule  and  the  capacity  rule  for 
k=l,2, . . . ,8  . 

In  the  printout,  ;iIND(J)  represents  defined  by  expression  (46), 

and  MINSD  represents  X  •  As  shown  in  the  printoit,  MINSD  is  0  for 

.1  ^ 

k=l,2,...,8  ,  and  so  t'ne  facility  usage  rule  does  not  'orce  any  facilities 

into  usage;  and  as  shown  by  matrix  CX(1,J)  for  k=l,2,.,.8  ,  the  capacity 

rule  does  not  fax  any  x.,  variables. 

■LJ 

Step  3;  Solve  the  relaxed  problem  (PR^)  . 

In  the  printout  FLB(K)  represents  ,  given  by  expression  (34), 

for  computing  FC.  ,  and  C(I,J)  represents  c..,  defined  by  expression 

IJ  X' 

(38).  Being  at  the  root  node,  X.  =  1.  Further  the  solution  of  problem 
(PRp)  ,  i.e.,  X  =  {x..}  is  shown  in  the  printout  by  SOLX(J)  which 

for  (PR.)  is  X  =  {x,  ,  =  x,T  =  x-.,  =  x,  ,  =1}  . 

1  41  42  2J  AA 

Step  4:  Compute  the  lower  bound. 

The  expressions  (40)  and  (41),  i.e., 

LOWS  =  zd’R^)  ^“^0) 

=  c...  +  X  min  c .  -f  FC  , 

.  ..  ij2  .  ^  (^1) 

jeW  jr.W  1 


(i, j)es 
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are  represented  in  the  printout  as 

LOWB  -  MINSC  +  FC 

=  729839.3125  +  0  =  729839.3125 


Step  5:  Compare  I.OWB  with  BUBS. 

Since  LOWB  <  BUBS,  go  to  Step  6 

Step  6:  Check  if  solution  X  is  feasible  in  problem  (P),  i.e., 
expression  (42)  is  satisfied. 


j ;  ‘‘ijt  »ij  i  Vk 

X.  .eX 
ij 


Vk 


(42) 


In  the  printout,  NSUKD  represents  the  left-hand  side  of  this  inequality, 
and  for  each  k  ,  the  capacity  constraints  are  satisfied. 

Step  7:  Compute  the  upper  bound. 

UPB  is  given  by  expression  (44),  i.e., 


UPB  =  X  E  a  X  +  Z  by  (44) 

i  j  k 

x.,CX 

ij 


In  the  printout,  the  corresponding  expression  is  represented  as 

UPB  =  NSUMA  +  FCUB 

=  678,502  +  101,000  =  779502.0  . 

Step  8:  Compare  UPB  with  BUB. 

Since  UPB  <  BUB,  go  to  Step  9. 

Step  9:  Set  BUB  =  779502.0  .  Since  EPS  =  0.0,  BUBS  =  BUB. 

Since  LOWB  <  BUBS,  go  to  Step  10. 

Step  10:  Left  branching  is  possible  since  W  =  <p  as  shot/n  by 
vector  FIX(J);  go  to  Step  11. 

Step  11:  Apply  the  bounding  rule  and  the  branching  rule. 


ij 


According  to  our  bounding  rule,  a  free  x 


variable  is  excluded 
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(fixed  at  0)  for  which 

(c.  1  -  c.  .„)  >  (BUBS  -  LOWB)  for  jeW  and  (i,j)eS  (51) 

IjX, 

For  ,  (210,381,4375-145,201.5)  >  (779, 502 . 0-729, 839. 3125) . 

This  also  holds  for  x^.  and  x, ,  ,  i.e.,  the  bounding  rule  results 

33  14 

in  fixing  ,  x^^  ,  and  x^^  at  0.  This  is  shown  in  the  printout 

by  matrix  CX(I,J)  where  the  corresponding  variables  have  been  assigned 
the  value  2  because  of  the  bounding  rule. 

The  branching  rule  directs  us  to  select  a  free  x..  variable 

°  ij 

corresponding  to  c.  .o  for  which  D,,  =c.  .n  -c.  .p  is  maximized 

over  all  j  .  In  the  printout,  c.  ,  c.  ,  and  D.„  are  represented 

by  NMINC(J),  MINC(J)  and  DIFBR(J),  respectively.  Since  D2J  is  the  maxi¬ 
mum,  x^2  is  selected  as  the  next  left  branching  variable.  This  is  shown 
in  the  printout  by  BRl  and  is  represented  as  (100  it-.')  e.g.,  402  . 


Using  the  terminology  employed  in  Figures  4a  and  4b,  the  x^^ 

variables  fixed  at  0  or  1  in  the  partial  solution  will  be  shown 

as  Sj^  =  {-  103,  -  303,  -  104,  402}  .  In  the  computer  printout,  vector 

STX  displays  the  x..  variables  fixed  at  0  or  1.  The  representation 

of  the  variables  is,  however,  somewhat  different.  An  x..  variable 

11 

fixed  at  0,  due  to  any  rule,  is  shown  as  -  (100  i  +  j)  -  1,000,000, 
e.g.,  x,T  is  shown  as  -  1,000,103;  an  x..  variable  .  ixed  at  1  due 

to  the  brandling  rule  is  represented  as  (100  i  +  j),  e.g,,  x^2  402; 

and  an  x..  variable  fixed  at  1  due  to  a  rule  other  than  the  branching 
ij 

rule  is  shown  as  (100  1+1)+  1,000,000,  e.g.,  X2,  is  represented  as 

1,000,203. 

In  the  printout,  vector  STX  represents  updated  partial  solution 


We  now  move  Co  Node  2. 


Node  2 

The  updated  matrix  CX(I,J)  and  vector  FIX(J)  are  displayed  in  the 
printout . 

Step  2:  Apply  the  facility  usage  rule  and  the  capacity  rule  for 
k=i,2,...,8  . 

As  shown  in  the  printout,  MINSD  (representing  ^  ’  being 

j  ^ 

positive  for  I =1,2, 3, A,  and  5  ,  these  facilities  are  forced  into 
usage.  Further,  for  k=A  ,  expression  (47)  holds  for  and  , 

i.e.. 


(180-0)  >  (200-30),  and 

(180-0)  >  (200-30),  respectively. 

As  shown  by  matrix  CX(I,J)  in  the  printout,  these  two  variables  are 
excluded  (fixed  at  0)  by  Che  capacity  rule.  Since  the  capacity  rule 
results  in  fixing  at  least  one  variable  in  the  first  cycle,  another 
cycle  is  repeated  as  displayed  in  the  printout.  The  second  cycle 
does  not  fix  any  more  variables.  Vector  STX  is  updated  accordingly. 

Step  3:  Solve  the  relaxed  problem  (PR2)  . 

6j^2  represented  by  FLB(K)  ,  c^j2  represented  by  matrix  C(I,J), 
and  solution  X  represented  by  SOLX(J)  are  displayed  in  the  printout. 

Step  4:  Compute  LOWB. 

LOWS,  from  the  printout,  is  equal  to  749011  4375. 

Step  5:  Compare  LOWB  with  BUBS. 

Since  LOWB  <  BUBS,  go  to  Step  6. 

Step  6:  Check  if  solution  X  is  feasible  in  (P) . 

In  the  printout,  for  k=4  ,  NSUI'ID  =  290  >  200  ,  i.e.,  expression 
(42)  is  not  satisfied,  and  we  go  to  Step  10. 

Step  10:  As  shown  by  vector  FIX(J),  left  branching  is  possible 
and  we  go  to  Step  11. 
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Step  11:  Apply  the  bounding  rule  and  the  branching  rule. 

As  displayed  by  matrix  CX(I,J)  in  the  printout,  the  bounding 
rule  results  In  fixing  X22^  and  at  0.  Now,  for  j=4  ,  except 

for  x^^  ,  all  the  x^^  variables  are  fixed  at  0;  therefore  x^^  , 

is  fixed  at  1.  This  is  reflected  by  matrix  CX(I.J),  and  vector  FIX(J). 
Vector  STX  is  updated  accordingly. 

The  branching  rule  selects  x^^^  as  the  next  branching  variable. 

This  is  shown  in  the  printout  by  BRl,  and  vector  STX  is  updated 
accordingly. 

We  now  n.3ve  to  Node  3. 


Node  3: 

The  updated  matrix  CX(I,J)  and  vector  FIX(J)  are  displayed  in  the 
printout. 

Step  2:  Apply  the  facility  usage  rule  and  the  capacity  rule  for 
k=l,2,...,8  . 


The  faciiity  usage  rule  forces  facilities  1  to  5,  and  8  into 
usage.  For  k=4  ,  the  capacity  rule  excludes  and  ,  i.e., 

fixes  them  at  0;  and  for  j=3  ,  all  but  being  fixed  at  0,  X22 

is  fixed  at  1.  This  is  displayed  in  the  printout  h>  matrix  CX(I,J) 
and  vector  FIX(J).  Vector  STX  is  updated  accordingly. 


Although  the  capacity  rule  has  fixed  at  least  one 


ij 


variable 


during  the  Initial  cycle,  another  cycle  is  not  necessary, 
by  vector  FIX(J)  which  represents  set  W,  since  we  have  an 


as  displayed 


ij 


variable 


fixed  at  1  for  each  of  the  n  columns  (activities) . 


Step  3:  Solve  the  relaxed  problem  (PR^)  . 

SOLX(J)  displays  the  solution  for  the  relaxed  problem. 


Step  h:  Compute  LOWS. 


LOWB,  shown  in  the  printout,  is  equal  to  779502. 0  . 
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Step  5:  Compare  LOWB  with  BUBS. 

Since  LOV/B  =  BUBS,  go  to  Step  12. 

Step  12:  Backtrack. 

We  backtrack  by  moving  to  the  parent  Node  2,  and  branching  to  the 
right  by  setting  =  0  (since  the  right  branch  has  not  yet  been  explored). 

In  the  printout,  this  is  accomplished  by  observing  the  last  entry  in  vector 
STX,  and  moving  backwards,  one  entry  at  a  time,  until  we  find  a  positive 
entry  without  1,000,000  added  to  it.  The  corresponding  x^^  variable 

is  fixed  at  0,  and  we  move  to  the  right  branch  node.  Matrix  CX(I,J),  vector 
FIX(J)  and  vector  STX  are  updated  accordingly.  As  disnlayed  in  the  print¬ 
out,  entry  401  in  vector  STX  is  such  an  entry,  and  variable  x^^^  is  fixed 

at  0  for  branching  to  the  right.  This  is  shown  in  the  printout  by  BRO  as  401. 
The  updated  vector  STX  is  also  displayed. 

We  now  move  to  Node  4. 


Node  4 

The  updated  matrix  CX(I,J)  and  vector  FIX(J)  art  displayed  in  the 
printout , 

Step  2:  Apply  the  facility  usage  rule  and  the  capacity  rule  for 

k=l,2,...,b  . 


As  displayed  in  the  printout,  for  k=4,  MINSr=230  >  200  ,  i.e., 
the  right-hand  side  of  inequality  (47),  (Sj^  -  Z  djj^j^)  <  0  ,  and 

according  to  our  backtracking  rules,  we  backtrack,  i.e.,  go  to  Step  12. 


Step  12;  Backtrack. 


We  backtrack  to  tlie  parent  Node  2,  and  since  the  right-hand  branch 
has  already  beep  explored,  backtrack  to  Node  1  and  to  the  right-hand 
branch  by  fixing  x^2  ®  •  This  is  shown  in  the  printout  by  BRO  as  402, 

and  vector  STX  is  updated  accordingly. 

We  now  move  to  the  next  node,  i.e..  Node  5. 


Branch-and-Bound  Tree 


We  continue  the  branch-and-bound  procedure  from  one  node  to 
another  until  we  backtrack,  to  the  root  node  and  find  that  the  right 
branch  has  already  been  explored.  The  procedure,  then,  terminates  and 
the  solution  corresponding  to  the  best  upper  bound  is  the  optimal  solu¬ 
tion. 

For  this  problem,  a  total  of  nine  nodes  are  explored  and  the 
optimal  value  equals  779502.0.  The  optimal  solution  is  ^42~ 

"  ^4  "  ^  =  y2  “  ^3  "  ^4  "  ^5  "  ^8  ^ 

displayed  in  the  computer  printout  on  the  last  page  of  Appendix  B. 

Figure  7a  presents  the  branch-and-bound  tree  ff r  this  problem, 
and  shows  the  node  numbers,  the  bounds,  and  the  bra.iching  variables. 

In  order  to  demonstrate  the  role  of  the  capacity  rule  and  the 
bounding  rule,  Figure  7b  displays  the  variables  fixed  as  0  or  1 

by  these  rules  for  this  test  problem. 

The  cumulative  effect  of  the  various  rules,  including  the  facility 

usage  rule,  the  capacity  rule,  and  the  bounding  rule,  .lakes  the  branch- 

and-bound  procedure  quite  efficient.  Further,  the  sto  age  and  updating 

of  the  X  variables  fixed  at  0  or  1  is  done  in  a  manner  that  makes 
ij 

utmost  use  of  ti.e  relevant  information  at  the  preceding  node. 
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Figure  7a 

Branch-and -bound  tree  for  a  test  problem 
(Test  Problem  with  m=5,  n=4,  and  p=8) 


Node 

Capacity  Rule 

Bounding  Rule 

1 

^13=°  ’  ^^33=^  •  "14=° 

2 

^34=0  .  ^54=° 

X2i=0  ,  X2^=0  ,  x^^=l 

3 

x^3~0  f  *  ^23”"^ 

5 

^3°®  ’  ^34=° 

6 

Xj^j=0  ,  X2^=0  ,  x^j^=0  ,  X22=0  1  ’'54^*^ 

8 

9 

>^53"’ >  ^44^'^’  ’^54"®’  ’^24"' 

Figure  7b 

Variables  fixed  by  the  capacity  rule  and  the  bounding  rule 


5.  COMPUTATIONAL  TEST  RESULTS 


The  computer  program  ZIPCAP  has  been  tested  on  several  problems. 
Although  primarily  designed  for  capacitated  problemr.  (i.e.,  where  the 
capacity  constraints  are  active),  the  program  can  also  be  used  for  solving 
uncapacitated  problems  as  a  special  case.  Since  the  dita  available  for 
capacitated  problems  were  limited,  some  uncapacitati-.d  problems  were  also 
considered  foj  testing  the  program.  (Most  of  the  data  were  furnished  by 
Professor  Pinkur  and  are  related  to  his  work  on  multi-echelon  Inventory 
systems. ) 

Table  4  presents  the  test  results  of  ZIPCAP.  In  order  to  verify 
the  optimal  solutions,  the  test  problems  were  also  solved  by  using  the 
0-1  Integer  programming  code  RIP30C  [Geoff rion  and  Nelson  (1968)], 

In  the  table,  the  problem  size  shows  the  number  of  designs  (ra) , 
activities  (n),  and  facilities  (p) .  This  is  equivalent  to  solving  a 
problem  having  ran+p  variables  and  n+p  constraints.  The  elapsed  time 
represents  the  time  in  seconds  to  solve  the  problem,  excluding  the  time 
to  read  and  write  the  input  and  to  write  the  outpu,:.  The  total  number 
of  nodes  explored  by  ZIPCAP  for  a  specified  set  of  options  is  also  shown. 

Both  RIP30C  and  ZIPCAP  were  run  on  an  IBM  3031  al  The  George 
Washington  University.  The  last  problem  in  the  table  was  not  run  using 
RIP30C  because  of  the  code's  capacity  limitation  to  90  variables  and 
50  constraints. 

The  test  problem  with  m=3  ,  n=4  ,  and  p=5  has  three  variations, 
using  different  values  for  the  facility  capacities.  The  data  for  the 
variable  costs  a^j  ,  fixed  cost  bj^  ,  and  the  capacity  requirements 

d...  are  given  in  the  other  document,  i.e.,  Chhabra  and  Soland  (1980). 

Ij  k 

For  the  test  problem  with  m*5  ,  n=4  ,  and  p^S  ,  runs  4a,  4b, 
and  4c  are  the  same  except  for  nc  different  intermediate  steps'  option 


RIP30C  ZIPCAP 


19.397  277 


(ISTEP)  and  this  results  in  slight  differences  in  the  time  taken  to  solve 
the  problem.  Runs  4a  and  4d  differ  in  that  4d  does  nJt  use  the  capacity 
rule;  the  resulting  difference  in  the  total  number  of  nodes  explored 
to  reach  the  optimal  value  points  to  the  effectiveness  of  the  capacity 
rule  in  conjun'tion  with  the  bounding  rule. 

Run  5  shows  the  results  for  an  uncapacitated  problem  with 
m=10  ,  n=8  ,  and  p=8  .  Option  ICAPR  is  not  usid  since  the  capacity 
rule  is  not  useful  for  an  uncapacitated  problem. 

Another  uncapacitated  problem  with  m=10  ,  n-30  ,  and  p=8  is 
solved  in  runs  6a  and  6b.  In  run  6a,  the  epsilon  value  (EPS)  is 
specified  as  0.002.  The  solution  value  found  by  exploring  125  nodes 
may  be  subopt  Imal  but  is  guaranteed  to  be  within  +0,2  percent  of  the 
optimal  solution  value.  Run  6b  is  made  with  an  epsilon  value  (EPS) 
of  0.0,  and  the  optimal  solution  value  is  found  in  277  nodes.  A  comparison 
of  runs  6a  and  6b  shows  that  the  number  of  nodes  ;,s  less  than  half  for  a 
solution  value  that  may  be  suboptimal  but  very  close  to  the  optimal  solu¬ 
tion  value  as  compared  to  the  number  of  nodes  for  an  optimal  solution 
value. 

In  general,  a  small  difference  between  a  solution  value  that  may 
be  suboptimal  ani  the  optimal  solution  value,  translates  into  a 
significant  difference  in  the  corresponding  number  of  nodes  and  the 
solution  time  required. 


6.  FURTHER  CONSIDERATIONS 


It  was  mentioned  in  Chapter  2  that  it  is  possible  to  consider 
alternative  formulations  of  problem  (P) ,  and  also  to  c  )nsider  choices  of 
Lagrange  multipliers  other  than  =  bj^  with  the  puroose  of  obtaining 

"tighter"  bounds  which,  in  turn,  would  further  improve  the  efficiency 
of  the  branch-and-bound  procedure.  These  aspects  will  be  discussed  in 
this  Chapter. 


6.1  Alternative  Formulations 

Problem  (P)  can  be  reformulated  by  adding  additional  constraints 
such  that  the  corresponding  Lagrangian  relaxation(s) ,  if  solved,  would 
provide  "tighter"  bounds.  If  such  a  relaxation  does  not  possess  the 
integrality  property,  then  it  provides  an  equal  or  better  bound  compared 
to  that  from  an  LP  relaxation,  as  mentioned  in  Chapter  2. 

Two  alternative  formulations  of  problem  (P),  along  with  their 
Lagrangian  relaxations,  are  given  below. 
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Since  =  1  or  0,  each  constraint  of  (52)  is  either  equivalent  to 

^  (if  e^j^  =  1)  or  else  is  redundant  (if  =  0).  Problem  (API) 

thus  has,  at  ntost,  mnp  additional  constraints  relative  to  problem  (P) . 

Two  Lagrangian  relaxations  are  now  considered  for  problem  (API). 

The  first  Lagrangian  relaxation  is  obtained  with  respect  to 
constraints  (5')  by  introducing  nonnegative  Lagrange  multipliers 
Vj^  ^  0  to  get 


Minimize 


E  Z  a  X..  +  I  b  y  -  5^  v  /y  -EE  r  .  x.  \ 
^  j  ij  ij  k-'k  ^  k  ^  ^  ijk  ij  j 

subject  to  (2),  (52),  and  (6),  or  equivalently. 


(ALRI^)  t 


Minimize 

EEx..  /a  +Ev  r.  ) 

i  j  ij  ^  ij  ^  k  ijkj 

-  1  \  (\  - \) 

(53) 

subject  to 

Ex.,  =1 
i 

Vi 

(2) 

"ik^jl  ^k 

Vi, j ,k 

(52) 

X..  ,  y  =  0  or  1 

Vi, j ,k 

(6) 

Another  Lagrangian  relaxation  of  problem  (API)  is  obtained  with 
respect  to  constraints  (5')  and  (52)  by  introducing  noni egative  Lagrange 
multipliers  v  and  i  respectively,  to  get 

K  IJ  K 


Minimize  E  E  a .  .  x .  .  +  E  b,  y, 

i  j  ij  ij  k  k 


-  Z  V, 


(>'k  -  “  'ijk  -ij) 


-  ^  “  *ijk  (^k  -  'ik  "ij) 


subject  to  (2)  and  (6),  or  equivalently. 
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(ALRl  ,)  X 

V,  A 


Minimize 


Z  Z  X.  ,  /a..+Zv  r.  +Ie  ^..,1 
i  j  AJ  \  k  k  ijk  ^  :.k  ijkj 


Subject  to 


k 

'ijk 

A  1  J 

(54) 

I  X  .  , 

=  1 

Vi 

(2) 

i 

Xij  , 

y^^  =  0  or  1 

Vi,  j  ,k 

(6) 

For  this  problem,  the  solution  is: 


and  X . . 
ij 


y  =  0 

if  (v 

+  I  I  X.  -  b,  \ 

1  0  , 

■'k 

[  k 

i  j  V 

=  1 

“  (a 

+  I  Z  X.  -  b  \ 

^  .  IJk  k) 

1  V 

o 

if  i  minimizes 

1  a.  .  +  Z  V,  to 

1  id  k 

k  'i" 

over  Si 


We  need  good  choices  of  Lagrange  multipliers  Vj^  with  which  t;o 
solve  problem  (ALRl^) ,  and  of  Lagrange  multipliers  Vj^  and 
which  to  solve  problem  (ALRl  ,).  Problem  (ALRl  )  d'^es  not  possess  the 

V  ,  A  V 

integrality  property,  thus  offering  the  hope  of  a  tigi.t  bound,  but  has 
more  constraints  and  is  difficult  to  solve  compared  to  problem  (ALRl^ 

which,  on  the  other  hand,  involves  more  Lagrange  mu  1 1  ip iers  . 


6.1,2  Alternative  Formulation  2 

Another  formulation  of  problem  (P)  is  similar  to  problem  (API) 
except  for  a  modification  in  constraints  (5'),  i.e.. 
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^  Minimize 


Subject  to 


(AP2) ^ 


V 


"ij  "ij  ^  1  w 

(4) 

X.  .  =  1 

IJ 

VJ 

(2) 

r . . ,  X, .  <  1 

ij  k  ij- 

Vk 

(55) 

e ,,  X .  .  <  y, 

ik  ij  —  ■'k 

Vi, j  ,k 

(52) 

X . .  ,  y,  =  0  or  1 
ij  -^k 

Vi,  j  ,k 

(6) 

A  Lagrangian  relaxation  with  respect  to  constrain  . s  (55)  and  (52)  is 
obtained  by  introducing  nonnegative  Lagrange  multipliers 
to  get 


and  X 


ijk 


Minimize  Z  Z  a  .  ,  x .  .  +  Z  b,  y, 

i  j  ij  xj  ^  k  k 


-  Z  V, 

k 


'k 

Subject  to  (2)  and  (6),  or  equivalently, 


(ALR2^^^)i 


Minimize 


Z  Z  X..  (a..  +  Z  y,  r  +  Z  e  X 

i  j  ij  ^  ij  k  k 


1 5  hjk) 


Subject  to 


Z  X ,  .  =  1 

,•  xj 


x^  ,  ,  y^^  -  0  or  1 


^  ''k 

k 

Vi,  J  ,k 


(56) 

(2) 

(6) 


For  this  problem,  the  solution  is: 


yfc  =  0  if  Z  Z  X^jj^  ^  bj^  , 


i  j 

z  >: 

i  j 


^  >■  \jk  i  \  • 


1  if 


69  - 


and  X, .  =  1  if  i  minimizes 
Ij 


(‘y  ^  "k  'y-  k  “ik  *yk - 


Here  again,  we  need  good  choices  of  the  Lagrange  multipliers 

and  A^.,  with  which  to  solve  problem  (ALR2 

ijk  v.A' 

6.1.3  Choice  of  Lagrange  Multipliers 

Each  of  the  relaxations  (ALRl  , )  and  (ALR2  , )  involves  p  v, 

v,A  v,A  k 

Lagrange  multipliers  and  mnp  A...  multipliers.  If  we  have  good  choices 

ij  k 

of  these  multipliers,  the  resulting  solutions  of  the  relaxed  problems 

should  provide  "tighter"  bounds  (because  of  the  additional  constraints) 

than  the  bound  from  relaxation  (LR  ).  Since  relaxa;ions  (ALRl  -, )  and 

V  v,A 

(ALR2  ,)  are  similar  to  a  great  extent,  only  the  relaxation  (ALRl. 

V, A  V  >  A 

will  be  considered  for  further  discussion. 

By  looking  at  expression  (54)  of  the  formulation  (ALRl  ^),  a 
meaningful  choice  of  the  Lagrange  multipliers  and  appears 

to  follow  from  setting 


V,  +  Z  Z  A. =  b,  Vk 

k  ^  j  ijk  k 

so  that  each  of  the  A,.,  can  be  chosen  as 

ijk 


^k  “  '"k 
^ijk  ~iifz  e,“  \ 


if  -  1  . 


0  otherwise 


The  solution  for  problem  (ALRl  ,)  is  then  to  select >  from  each  column  j 

V,  A 


an  x^^  variable  which  minimizes 


over  I 


Arbitrary  values  were  considered  for  the  v^^  (e.g.,  Vj^  equal  to 

3/4  bj^  ,  1/2  bj^  ,  1/4  bj^  ,  and  0)  ,  the  were  then  computed  from 


}  •  •*  xV*. 
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(58),  and  the  test  problem  with  three  designs  (m) ,  four  activities  (n) 
and  five  facilities  (p)  was  solved.  Three  cases  with  different  capacities 
Sj^  (as  specified  in  Chapter  5,  Table  4)  were  tried  for  the  solutions  at 

the  initial  node.  The  results,  however,  were  not  cone '.usive  in  terms  of 
providing  a  meaningful  choice  of  the  Lagrange  multipliers  Vj^  (and  of 


Since  the  relaxation  (ALRl  possesses  the  integrality  property, 

V,  A 

a  choice  of  the  multipliers  as  the  optimal  values  of  the  dual  variables 
of  the  corresponding  linear  programming  problem  woulo  provide  a  solution 
as  good  as  the  LP  solution  (as  stated  in  Chapter  2) .  Uc  do  not  propose  to 
solve  linear  programs  as  a  part  of  our  branch-and -bound  methodology. 
However,  we  have  made  some  LP  runs,  basically  to  see  if  the  results  provide 
insight  leading  co  the  choice  of  the  Lagrange  multipliers,  and  also  to 
see  if  the  resulting  LP  solutions  are  "close"  to  the  integer  solutions. 
These  results  are  given  below. 


The  LP  formulation  (API)  corresponding  to  prohiem  (API)  is: 


Minimize 


Subject  to 


V 

The  constraints  x.. 

ij 


I  E 

i  j 

I 

i 

E  E 


a, .  X. .  - 
id  ij 

k 

(A) 

X.  .  =  1 
ij 

Vi 

(2) 

ij 

Vk 

(5’) 

^ik  '‘ij 

Vij.k 

(52) 

<  1 

Vk 

(15) 

>  0 

Vi, j ,k 

(16) 

are  implicit  in  constraints  (2), 


p=5 


Problem  (API)  was  solved 
and  three  different  cases 


for  the  test  problem  w.. ch 
for  the  capacities  s^^  (as 


m=3 ,  n=4 , 

specified 


and 

in 
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Chapter  5,  Table  4).  Each  case  was  solved  using  the  IMSL  (International 
Mathematical  and  otatistical  Library)  Code  ZX3LP  on  an  IBM  3031  at  The 
George  Washington  University. 

Note  that  the  formulation  (API)  has  up  to  mnp  more  constraints  than 
the  LP  formulation  (P)  given  in  Chapter  2.  For  our  test  problem,  this 
translates  into  solving  a  problem  of  17  variables  and  50  constraints 
corresponding  to  formulation  (API)  as  against  17  variables  and  14  con¬ 
straints  corresponding  to  formulation  (P) . 

Table  5  lists  the  solution  values  for  each  of  the  three  cases  with 
different  capacities  for  the  small  problem  with  thr^^e  designs,  four  activi¬ 
ties  and  five  facilities.  The  solutions  to  problem!  (P)  and  (API), 
obtained  from  ZX3LP,  show  the  optimal  solution  values,  the  optimal  values 
of  the  variables  and  y^^  ,  and  the  optimal  values  cf  the  dual 

variables  corresponding  to  the  Lagrangian  relaxations  (LR  )  and  (ALRl  . ) , 

V  V ,  A 

l.e.,  Vj^  associated  with  the  capacity  constraints  (j')  and 

associated  with  the  constraints  (52).  The  table  also  shows  Z(P),  and 
the  Lagrangian  solution  value  Z(LR^)  obtained  by  setting  Vj^  =  bj^  for 

all  k  at  the  root  node,  i.e.,  Z(LR,). 

b 


As  expected,  the  LP  solutions  for  each  of  the  three  cases  show 
Z(APl)  to  be  considerably  higher  than  Z(P),  and  closer  ts  Z(P),  thereby 
providing  a  tighter  bound.  As  for  the  Lagrange  multipliers  Vj^  and 


X 


ijk 


the  following  relationships  are  observed. 


z  z 

i  j 

-1, 

1  \ 

V  k  ,  and 
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i  j 

^ijk  -  ^k 

Vk  . 

Also , 
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< 
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°  ^  ?  ^ijk  ° 

1  j 

b,  ,  and 
k 
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'’k  •  “  *ijk  - 
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Although  the  relationship  among  various  X,  values  is  not  apparent, 

ij 

the  above  observations  are  useful  in  further  exploring  some  good  choices 

of  the  Lagrange  multipliers  for  the  relaxation  (AI,R1  ,),  as  discussed 

V ,  A 

earlier. 

As  for  the  "closeness"  of  the  LP  solution  to  that  of  the  integer 

solution,  most  of  the  solution  values  x  and  y  of  problem  (API) 

ij  k. 

fractional,  and  their  rounding  off  to  0  or  1  does  not,  in  general, 
seem  to  correspond  to  the  optimal  integer  solution  values  .  and 

of  problem  (P). 

Table  5  also  displays  Z(LRj^)  at  the  root  node  for  each  of  the 
three  cases.  For  Sj^  =  3000  ,  Z(LRj^)  =  Z(P)  ,  and  the  Lagrange  multi¬ 
pliers,  as  reflected  by  the  values  of  the  dual  variables  of  problem  (P) , 
are  equal  to  bj^  for  all  k  .  This  is  expected  from  Theorem  3  and  our 

discussion  of  the  integrality  property  in  Chapter  2.  Further,  the  dual 
variables  of  (P)  for  the  first  two  cases  (i.e.,  Sj^  =  7C0  Vk  ,  and 

s  =  400,...)  have  values  v,  >  b,  from  Theorem  1. 

The  Z(LRj^)  values  in  Table  5,  however,  take  no  consideration  of 
the  capacity  rule  and/or  the  facility  usage  rule  of  our  branch-and- 
bound  procedure.  These  rules,  by  fixing  certain  x^^  values  to  0 

or  1,  and  by  forcing  certain  facilities  into  the  solution,  provide 

an  improved  lower  bound.  As  per  our  branch-and -bound  procedure  the 

improved  lower  bound  at  the  root  node  is  obtained  by  solving  problem 

(PR  ) .  For  example,  for  s,  =  700  Vk,  the  values  of  Z(LR,)  and  Z(PR, ) 

X  K  b  1 

are  shown  in  Figure  8.  The  figure  also  shows  the  values  of  Z(P) , 

Z(APl),  and  Z(P).  The  branch-and-bound  procedure  rules  provide  an 

improved  value  of  the  lower  bound  Z(PRj)  compared  to  Z(P).  It  appears 

that  some  good  values  of  the  Lagrange  multipliers  of  the  relaxation 

(ALRl  ,),  if  found,  could,  in  conjunction  with  these  rules,  provide 
V ,  A 

significant  improvement  over  Z(APl),  and  without  the  need  to  solve  an 
LP  problem. 
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6.2  Subgradient  Method 

It  was  mentioned  in  Chapter  2  that  setting  the  Lagrange  multipliers 
equal  to  for  all  k  provides  a  good  starting  point  in  solving 

the  Lagrangian  relaxation  (LR^)  of  problem  (P) .  From  Theorem  3,  this 

choice  is  optimal  (in  terms  of  providing  the  tightest  lower  bound)  if 
the  resulting  solution  is  feasible  in  problem  (P) .  In  other  cases, 
i.e. ,  where  the  resulting  solution  is  not  feasible  Ir  problem  (P) ,  this 
choice  is  generally  not  optimal  and  it  is  possible  to  tighten  the  bounds 
by  considering  values  of  Vj^  ^  bj^  (from  Theorem  1)  .  Oi.e  method  that 

seems  useful  in  providing  such  a  choice  is  the  subgradient  method. 

The  subgradient  method  is  an  adaptation  of  the  gradient  method 
in  which  gradients  are  replaced  by  subgradients.  Through  a  heuristic 
choice  of  the  step-size,  this  method  has  been  successfully  used  to 
provide  Improved  bounds  and  sometimes  optimal  solutions  [for  details 
see  Held,  Wolfe,  and  Crowder  (1974),  Fisher  (1978),  and  Christofides 
(1980)].  The  fundamental  theoretical  result  is  that 

Z(LR  ^)^Z(D)  if  t*~^0  and  I  t9-^oo  as  g  — ♦  ®  , 

^  q=0 

where  t®  is  the  positive  step-size  t  for  the  g'.h.  iteration, 
and  Z(LR^®)  is  the  solution  value  of  the  relaxed  problem  (LR^)  using 

Vj^  values  obtained  at  the  gth  iteration. 

e+1 

In  the  case  of  problem  (P) ,  the  step-size  t®  ) or  iteration 
g  +  1,  given  that  we  have  a  solution  of  (LR^^),  is  given  by 


.8+1 


-  Z(LR  ®)] 


■  ^  J  '"ijk  ’‘ij 


8|  |2 


(59) 


2^1 

where  A  is  a  scalar  and  generally  between  0  and  2 ,  and 

Z  is  an  upper  bound  on  Z(LR^®),  frequently  obtained  by  applying 


a  heuristic  to  solve  problem  (P). 
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Given  the  vector  ot  multipliers 


,8+i 


is  generated  by 


Z  r 
i  j 


ijk 


(60) 


where  we  enforce 


in  our  case  of  problem  (P)  (because  of 


Theorem  1) . 


Justification  for  these  rules  and  conputationai  results  of 
applications  of  the  subgradient  method  are  given  in  Held  (1974). 

The  scalar  A  is  generally  initiated  by  setting  A^  =  2  and  halving 
subsequent  A's  whenever  the  resulting  solution  value  has  failed  to 
increase  in  some  fixed  number  of  iterations.  This  rule  has  performed 
well  empirically  [Held  (1974)  and  Fisher  (1978)]. 

For  the  test  problem  with  three  designs,  four  activities,  five 
facilities,  and  the  capacities  s^  =  400,  S2  =  400,  “  1000, 

s^  =  400,  s^  =  400,  the  Lagrangian  solution  obtained  it  the  root  node 

by  setting  =  b^^  for  all  k  ,  i.e.,  the  solution  to  problem  (PR^) 

is  infeasible  it  problem  (P) ,  i.e.,  it  violates  the  capacity  constraints. 
It  seems  that  the  subgradient  method  could  be  uselul  in  considering 
V  ^  b  with  the  ultimate  purpose  of  obtaining  a  tighter  lower  bound, 

and,  depending  on  the  revised  sol ut jon(s) ,  possible  amprovement  in  the 
best  upper  bound.  Another  posslblity  could  be  to  first  arbitrarily 
increase  the  relevant  Vj^  values  by  a  small  percentage  of  the  bj^ 

values  and  then  solve  problem  (LR^) ,  liopefully  to  improve  the  lower 

bound;  and  thereafter  to  use  the  subgradient  methoil  for  obtaining 
subsequent  values  of  v^^  ,  and  tiglitening  the  bounds. 


Both  of  tlie  areas  discussed  above,  i.e.,  the.  consideration  of 
alternative  formulations  of  problem  (P),  and  the  appl'.eation  of  the 
subgradient  metliod,  and  their  combination,  seem  useful  for  continued 
researcl)  in  term;-,  of  further  improving  the  branch-ani-bound  procedure 
for  solving  r’u.  multiactivity  multifacility  capacity-constrained  0-1 
assignment  problems. 
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fortran  IV  G  LEVEL  21 
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21PCAP,  A  ZERO-ONE  INTEGER  PROGRAM  IS  DESIGNED 
TO  SOLVE  MULTIACTI VITY  MULT  I F AC  I L I TY  CAPACITY- 
CONSTRAINED  PROBLEMS  HAVING  VARIABLE  AND  FIXED 
COSTS.  IT  ALSO  SOLVES  UNC AI  AC  I TATED  PROBLEMS  AS  A 

special  case 

INTEGER  0(35,35,301,  A(35,35),  CX{35,3f),  E ( 35 , 30 ) , 

1  8(301,  BS0LX(35),  BS0LY{30),  FL8 ( 30  I , F I  X ( 35 ) , F I X I ( 3 5  I , 

2  FUa(30i,  S(30t,  S0LX(35),  STX(l22ri 

REAL  MINC(351,  NMINC(35) 

DIMENSION  C(35,35l,  0IFBR(35>,  KT2(35k,  HIND(3'>I 

-INTEGER  BRO,  BRl,  FC,  FCUB,  P 

real  LOWB,  MAXOIF,  HINSC 

••••»*»*««*««*0PTIQNS  available:  IINPT,  ICAPR,  ISTEP,  IUNCAP,£PS 
11NPT=1  IF  INPUT  LISTING  DESIRED;  0  OTHERWISE 
ICAPR=1  IF  CAPACITY  RULE  TO  6E  USED;  0  OTHERWISE 
1STEP»0  IF  LISTING  OF  INTERMEDIATE  STEPS 
NOT  DESIRED.  IST£P=1  IF  SUMMARY  OF  BRANCH  £ 

BOUND  NODES  DESIRED.  ISTEP=2  IF  DETAILED 
LISTING  OF  intermediate  STE^S  DESIRED. 

1UNCAP  =  1  IF  SOLVING  AN  UNC AI' AC  I TATED  PROBLEM, 

0  OTHERWISE. 

EPS=  A  fractional  VALUE  IF  SUBOPTIMAL 
SOLUTION  DESIRED,  E.G.,  EPSILON  AS  0.005 
IMPLIES  SOLUTION  TO  BE  WITHIN  vC.5  PERCENT 
OF  THE  OPTIMAL  SOLUTION.  £PS=0.0  IF  OPTIMAL 
SOLUTION  DESIRED. 

£T=  ELAPSED  TIME  IN  SECONDS,  IF  SOECIFIEO,  AT 
WHICH  THE  NODE  AND  BOUNDS  RELATED  INFORMATION 
IS  PRINTED.  This  is  useful  IN  A  SITUATION  IF 
ISTEP=0  AND  THE  PROGRAM  TERMINATES  BEFORE 

reaching  the  final  solution. 

»««•> ad  input  DATA*********»***aa**««**»*********** 

read  10,  IINPT,  ICAPR,  ISTEP,  lUNCAP,  EPS,  ET 
10  FORMAT  (411,  F6.5,  F10.31 

M*  NUMBER  OF  DESIGNS 
N=  NUMBER  OF  ACTIVITIES 
P=  NUMBER  OF  FACILITIES 

READ  20,M,N,P 
20  FORMAT  (3151 

A(I,J|;  VARIABLE  COST  MATRIX 
READ  30,  ((A(1,JI,  I=1,M),J=1,NI 

30  FORMAT  (8I10I 

B(Kl:  FIXED  COST  V' CTOR 
READ  30,  (B(M ,K=1 ,P> 

IF  (lUNCAP.EQ.n  GO  TO  40 

S(KI:  CAPACITY  LIMIT  VECTOR;  REQUIRED  ONLY 
IF  IUNCAP=0 
READ  -.0,  (S(KI,K=1,P) 

0(I,J,Kl:  CAPACITY  USAGE  MATRIX;  REQUIRED 
ONLY  IF  IUNCAP=0 

DO  32  K=1,P 

READ  30, ( ( D( t,J,Kl ,  !>! ,M 1 , 1 , N 1 
32  CONTINUE 
DO  37  X=1,P 
DO  37  1  =  1, M 

IF  (0(  1 , 1 ,X I .EO.OI  GO  TO  35 
E(  I,KI  =  1 
GO  TO  37 
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00000  no 
00000120 
00000130 
00000140 
00000150 
00000160 
00000170 
00000180 
00000190 
00000200 
00000210 
00000220 
00000230 
00000240 
00000250 
00000253 
00000256 
00000260 
00000263 
00000266 
00000270 
00000273 
00000280 
00000290 
00000300 
00000310 
00000320 
00000330 
00000340 
00000350 
00000360 
00000370 
00000380 
00000390 
00000400 
000004 10 
00000420 
00000430 
0000Ch40 
00000450 
00000460 
00000470 
000004S0 
0000049C 
O0C0C50O 
00000510 
00000520 
00000530 
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0023 

35 

E(  I,KI-0 

000005*0 

002* 

37 

CONTINUE 

00000550 

» 

0025 

CO  TO  90 

00000560 

C 

E(I,KI:  DESIGN-FACILITY  MATRIX)  REOUIREO  ONLY 

00000570 

C 

IF  IUNCAP=1 

00000580 

c 

0026 

*0 

READ  *5,((EII,KI,I=1,M),K=1,P» 

00000590 

0027 

*5 

FORMAT  (80111 

00000600 

0028 

DO  80  K-1,P 

00000610 

o 

0029 

5  ( K  )  =N 

00000620 

0030 

DO  75  1=1, M 

00000630 

0031 

IF  (E(I,K).EQ.1I  GO  TO  65 

000006*0 

e 

0032 

00  60  J=1,N 

00000650 

0033 

0( I, J,K)-0 

00000660 

003* 

60 

CONTINUE 

00000670 

o 

0035 

60  TO  75 

00000680 

0036 

65 

00  70  J=1,N 

00000690 

0037 

0(  I,J,K)=1 

00000700 

o 

0038 

70 

CONTINUE 

00000710 

0039 

75 

CONTINUE 

00000720 

00*0 

80 

CONTINUE 

00000730 

0 

C 

000007*0 

00*1 

90 

PRINT  95,  IINPT,  ICAPR,  ISTEP,  lUNCAP,  EPS,  ET 

00000750 

00*2 

95 

FORMAT  (  •!•,'  OPTIONS  SELECTED  :  IINPT=',11, 

00000760 

o 

1 

L  '  ICAPR=',11,  •  ISTEP=',I1,  '  1UNCAP=',I1, 

00000770 

2  '  EPS=',F8.5,  •  ET=',  F10.3///) 

00000780 

00*3 

IF  ( IINPT. EQ. 01  GO  TO  168 

00000790 

c 

00** 

PRINT  100,M,N,P 

00000600 

00*5 

100 

FORMAT  ('O',  T55,  'INPUT  DATA', /IX,  T55,  ,//////lX , 

00000810 

1T*1«  'NUMbbR  OF  DESIGNS  *X,1*//1X»T*1 » 

2*NUMBER  OF  ACTIVITIES  (Nl«*,  IX,I*//IX,  T<.  I , 

3»NUNBER  OF  FACILITIES  <F>=*flX,  !*///> 

PRINT  105 

lOS  FORMAT  (  *Xt  'VARIABLE  COST  MATRIX  AIIiJ)'t/4Xt 
!■ 

DO  110  1^1, M 

110  PRINT  115,  I,  (A(I,Ji,J=l,N) 

115  FORMAT  ('O',  T6,  13,  *X,8I13,  *(/,  l*r,8I13)) 

PRINT  120 

120  FORMATC 'O' ,//*X, 'FIXEO  COST  VECTOR  B(K)',/*X, 

1'  .1— t,  /» 

PRINT  122,  (B(K),Kcl,P) 

122  FORMAT  ('O',  T15,  8113,  31/,  14X,8113»» 

PRINT  125 

125  F0RMAT('0',//*X, 'CAPACITY  LIMIT  VECTOR  S(I<I',/4X, 

I •  .  — »  ,  /  J 

PRINT  128,  (S(KI,Kxl,P) 

128  FORMAT  ('O',  T15,  3113,  3(/,  1*X,8I13)) 

PRINT  130 

130  FORMAT( 'O' ,//*X,  'CAPACITY  USAGE  MATRIX  0 ( I , J ,K > ' , /*X, 

1'-  - - ',/» 

00  150  Ksl,P 

PRINT  135, K 

135  FORMAT  ( '0 ' ,//5X , 'K  =  ' , 1 3/ » 

00  1*5  1=1, M 

PRINT  1*0,1  ,(0( I,J,KI ,  J  =  l,Ni 
1*0  FORMAT  ('O',  T6,  '1=',  13,  *X,8I13,  *(/,  1*X, 81131) 

1*5  CONTINUE 
150  CONTINUE 


00000820 

00000830 

000008*0 

00000850 

00000860 

00000870 

00000880 

00000890 

00000900 

00000910 

00000920 

00000930 

000009*0 

00000950 

00000960 

00000970 

00000980 

00000990 

00001000 

00001010 

00001020 

00001030 

000010*0 

00001050 

00001060 

OC001070 

00001030 

00001090 

00001100 

00001110 
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0069 

PRINT  155 

00001 120 

0070 

155 

FORMATi ‘O* i//4X, 'OeSIGN-f ACILITY  MATRIX  eiIfK!*,/4X, 

00001130 

( 

1»-.  ■■■■  '■■■■' 

00001140 

0071 

00  160  l>ltM 

00001150 

0072 

PRINT  158,  I,  (El  I  ,KI,K>1.P1 

00001160 

( 

0073 

158 

FORMAT  ('O',  T6,  •!=•,  13,  4X,8I13,  31/,  14X,S113n 

00001170 

0079 

loO 

CON . INUE 

00001180 

0075 

168 

IF  USTEP.EO.OI  GO  TO  190 

00001190 

:) 

0076 

IF  (ISTEP.EO.il  GO  TO  175 

00001200 

0077 

PRINT  170 

00001210 

0078 

170 

FORMAT  (  •0*,///55X, 'DETAILED  LISTING  OF  ITEPS',/) 

00001220 

,  \ 

0079 

GO  TO  190 

00001230 

0080 

175 

PRINT  180 

00001240 

0081 

180 

FORMAT  ( 'O', ///55X, 'SUMMARY  OF  STEPS',/) 

00001250 

o 

C 

00001260 

0082 

190 

auB*9999999. 

00001270 

0083 

8ues=  8ua/  (i.o«EPSi 

00001280 

o 

0084 

NSX^O 

00001290 

0085 

NOO=1 

00001310 

0086 

IBN0D=1 

00001315 

0 

0087 

INET=0 

00001320 

0088 

INSET=0 

00001330 

0089 

DO  205  J>:I,N 

00001390 

c 

0090 

f: X( J)=0 

00001400 

0091 

XT2( J)=0 

00001410 

0092 

DO  205  Ial,M 

00001420 

c 

0093 

CX(1,J)»0 

00001430 

0094 

205 

CONTINUE 

00001433 

0095 

LQ1=0 

00001436 

i  < 

0096 

LQ2=0 

00001440 

0097 

LR2=0 

00001443 

0098 

CALL  TIHETIITOJ 

000014„5 

0 

0099 

IF  (ISTEP.EO.OI  GO  TO  208 

00001448 

0100 

PRINT  220, NOO 

00001450 

0101 

208 

IFINSX.EO.OI  GO  TO  283 

00001453 

0 

C 

CX(I,J)  CONTAINS  FIXED  AND  FREE  XII, J)  VARIABLES. 

00001456 

c 

STXIINSI  CONTAINS  FIXED  XII, Jl  VARIABLES. 

00001460 

c 

CX(I,J)  AND  STXIINSI  ARE  UPDATED  BY  THE  CAPACITY 

00001480 

o 

c 

RULE,  THE  BOUNDING  RULE,  AND  THE  RULE  FOR 

00001490 

c 

BRANCHING  AND  BACKTRACKING. 

00001500 

c 

IN  CX(I,J)  A  FIXED  variable  .S  RECORDED  AS  1  OR 

00001505 

0 

c 

2,  AND  A  FREE  VARIABLE  AS  0. 

00001510 

c 

A  VALUE  OF  1  IMPLIES  THAT  THAT  PARTICULAR  VARIABLE 

00001515 

c 

IS  FIXED,  AND  FIXIJI  IS  SET  EOL'AL  TO  1  IMPLYING 

00001520 

0 

c 

THAT  COLUMN  J  HAS  A  FIXED  VARIABLE  OF  VALUE  1. 

00001525 

c 

A  VALUE  OF  2  IMPLIES  THAT  THAT  PARTICULAR  VARIABLE 

00001530 

c 

SHOULD  NOT  BE  CONSIDERED  FOR  CURRENT  COMPUTATIONS. 

00001535 

0 

c 

AN  X(I,J)  RECORDED  IN  CX(I,JI  AS  1  DUE  TO  THE 

00001540 

c 

BRANCHING  RULE  IS  RECORDED  IN  STXI.'NSI  AS  X*1004J. 

00001545 

c 

AN  XII, Jl  RECORDED  IN  CXII,J)  AS  1  DUE  TO  THE 

00001550 

(. 

c 

capacity  rule  OR  THE  BOUNDING  RULE  IS  RECORDED  IN 

00001555 

c 

STXIINSI  AS  (X*lOO*J)*100COOO. 

00001560 

c 

AN  X(I,J)  RECORDED  IN  CXII,J)  AS  2  IS  RECORDED  IN 

00001565 

c 

STXIINSI  AS  -IX*IOO*J I-IOOOOOU. 

00001570 

0102 

210 

IF  (IS^EP.tQ.O)  GO  TO  225 

00001580 

0103 

215 

PRINT  220, NOD 

00001590 

0104 

220 

format  I 'O', //6X, 'NODE  NUMBER',  17/1 

00001600 

C 

00001610 
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C 

BRO  IS  THE  RIGHT  BRANCHING  VARIABLE 

00001619 

0109 

229 

LX>BR0 

00001620 

0106 

IX>LX/100 

00001630 

0107 

JX*LX-IX*100 

00001640 

0108 

CX( lXtJX)-2 

00001690 

0109 

KT2UXI«KT2(  JX)4l 

00001660 

0110 

FIX! JXl^O 

00001720 

0111 

L01>LQ1-1 

00001729 

0112 

IF  (KT2IJXI.LT.(N-ltl  GO  TO  270 

00001730 

0113 

DO  299  1^1 fN 

00001740 

0114 

IF  (CXI  l.JX).EQ.2l  GO  TO  299 

00001790 

0119 

CX(I,JXI-1 

00001760 

0116 

NSXsNSX^l 

00001763 

0117 

STXINSXIa  ( I* 1004 3X141000000 

00001766 

0118 

FIX! JXI>1 

00001770 

0119 

LQl>LQl4l 

00001760 

0120 

FIXI IJXI«1 

00001790 

0121 

GO  TO  270 

0C001800 

0122 

299 

CONTINUE 

00001810 

0123 

270 

LQ2sO 

00001820 

0124 

LR2>0 

00001829 

0129 

GO  TO  283 

00001830 

0126 

272 

IF  (ISTEP.EC.OI  GO  TO  276 

00001840 

0127 

PRINT  220fNOO 

00001850 

C 

00001853 

C 

BRl  IS  THE  LEFT  BRANCHING  VARIABLE 

00001896 

0128 

276 

LQ2-0 

00001860 

0129 

LR2>0 

00001866 

0130 

LX«BR1 

00001870 

0131 

IX>LX/100 

00001879 

0132 

JX»LX-IX*100 

00001880 

0133 

CX(IX,JX)*1 

00001889 

0134 

FIX(JXI>1 

00001890 

0139 

LQ1bL0141 

00001892 

0136 

00  279  1=1, N 

00001699 

0137 

IF  ( lX.EQ.l)  GO  TO  281 

00001897 

0133 

279 

CONTINUE 

00001900 

0139 

261 

FIXll JXI=IX 

00001902 

0140 

283 

IF  (1STEP.NE.2I  CO  TO  303 

00001909 

0141 

289 

DO  299  1=1 ,N 

00001910 

0142 

PRINT  290,  I,(CX(I,JI,J=1,NI 

00001920 

0143 

290 

FORMAT  l/9X,<CXI I,JI*,4X,'1>*,13, 

2X,  2014/23X,  20I4I 

00001930 

0144 

299 

CONTINUE 

00001940 

0149 

PRINT  297, (FIXI JI,3a|,NI 

00001990 

0146 

297 

FORMAT  (/9X,'FIX(JI',12X,  2014/23X, 

20141 

00001960 

c 

**«*«*«****«««APP|.V  CAPACITY  RULE*** 

00001970 

c 

AND  UPDATE  CXII,JI  AND 

STXI INSI. 

00001980 

0147 

303 

DO  307  K=1,P 

00002000 

0148 

FLBIKI=0 

00002019 

0149 

307 

CONTINUE 

00002029 

0190 

310 

00  2000  K=1,P 

00002030 

c 

FIND  THE  SUM  OF  MINIMUM  DI1,J.X>  OVER 

EACH  J  FOR  A 

000U2040 

c 

GIVEN  K,  I.E.,  MINSC* 

SUM  OF  *IN0IJ1 

00002090 

0191 

MINSD-0 

OOOO2O0O 

0192 

DO  900  J=1,N 

00002070 

0193 

IFIFlXl JI.EO.OI  GO  TO  390 

00002030 

c 

IF  FIX(JI=1,  SET  MINU(JI=0(I,J,KI  FOR 

CX( t,Jl=l 

00002090 

c 

AND  MOVE  TO  NEXT  COLUMN  J 

00002 100 
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0154 

INOI^FIXK  J) 

00002110 

0155 

MINOlJl^OdNOltJ.Kt 

00002120 

( 

0156 

GO 

TO  800 

00002130 

0157 

350 

LK^O 

00002160 

0158 

1*1 

00002170 

c. 

0159 

MlNOl J)*0(I,J,KI 

00002180 

C 

SKIP  0(1. J.K)  WHEN  CX(I>J>*2  C  MOVE  TO  NEXT  RON  I 

00002190 

OlbO 

400 

1F(CX(I,J|.E0.2I  GO  TO  600 

00002200 

0161 

500 

’F(OIItJiK).LT.HINO(jn  HINDU  1*0 1  I.J.X) 

00002210 

01o2 

GO  TO  700 

00002220 

0163 

600 

LK*LK«1 

00002230 

0164 

IFd.GT.LKI  GO  TO  700 

00002240 

0165 

1*1*1 

00002250 

0106 

MINO(  JI*OdtJ,KI 

00002260 

o 

0167 

GO  TO  750 

00002270 

OIoS 

700 

1*1*1 

00002280 

0169 

750 

IFd.LE.MI  GO  TO  400 

00002290 

c 

0170 

800 

N1NS0*H1NS0*H1N0( J) 

00002300 

0171 

900 

CONTINUE 

00002310 

0172 

910 

IF 

dSTEP.NE.2)  GO  TO  960 

00002320 

o 

0173 

PRINT  950.  K.  MINSOt(MlNOUltJ*l,N) 

00002330 

0174 

950 

format  (•0*.'K,MINSD«(MINDUI.J*l,NI*.10Il0.4(/,44X.ail0) ) 

00002340 

0175 

960 

IF  (MINSO.EC.O)  GO  TO  975 

00002342 

0 

0176 

965 

IF(  FLBlKt.EO.ll  GO  TO  975 

00002344 

0177 

970 

F..B(KI»l 

00002346 

0178 

975 

IF  (lUNCAP.EO.l)  GO  TO  2000 

00002348 

o 

0179 

978 

IF  (ICAPR.EO.O)  GO  TO  2000 

00002349 

c 

FIND  BALANCE  AVAILABLE  CAPACITY  IBALO  FOR  A  GIVEN  K 

00002350 

c 

IF  IBALO  IS  NEGATIVE.  THEN  BACKTRACK. 

00002360 

L 

0130 

980 

IBALO-SIKI-MINSO 

00002380 

0181 

1000 

IF 

(IBALO. LT. 01  GO  TO  6200 

00002390 

0182 

00 

1500  J*l.N 

00002400 

0 

c 

SKIP  COLUMN  J  IF  F1X(3I*1 

00002410 

0183 

IF  (FIXUl.EQ.ll  GO  TO  1500 

00002420 

0184 

00 

1300  1*1, M 

00002430 

0 

c 

SKIP  RON  1  IF  CXd,JI*2 

00002440 

0135 

1100 

1F(CX( ItUI.E0.2l  GO  TO  1300 

00002450 

c 

COMPUTE  OIFFERENCE  BETWEEN  Dd.J.K)  AND  MIND(J). 

00002470 

o 

c 

IF  IT  IS  MORE  THAN  AVAILABE  BALANCE.  SET  CXd.JI*? 

00002480 

0186 

1200 

1D1F0*0(I.  J,K)>MINOU) 

00002490 

0197 

IF  (d01F0-lBAL0).LE.0>  GO  TO  1300 

00002510 

o 

0188 

CX'I,J)*2 

00002520 

0139 

NSX*NSX*1 

00002523 

0190 

STX(NSX )=-(!• 100*31-1000000 

00002526 

c 

LQ2  COUNTS  THE  NUMBER  OF  CXd.J)  VALUES  SET  EQUAL 

00002530 

c 

TO  2  IN  A  CYCLE 

00002540 

0191 

LQ2*Lg2*l 

00002550 

o 

c 

KT2(J)  KEEPS  AN  ACCOUNT  OF  CXd.JI  VALUES  SET  EQUAL 

00002560 

c 

TO  2  FOR  COLUMN  J 

00002570 

0192 

KT21 J»=KT2(J>*1 

00002580 

f 

c 

FOR  COLUMN  J.  IF  ALL  BUT  ONE  CXd.UI  VALUES  ARE 

00002590 

c 

EQUAL  TO  2.  SET  THAT  CXd.JI*!  C  SET  F1X(JI*1 

00002600 

0193 

IF(KT2( JI.LT.(M-1M  GO  TO  1300 

00002610 

1 

0199 

00 

1250  LR=1,M 

00002620 

0195 

IFICX(LR.JI.E0.2I  GO  TO  1250 

00002630 

0196 

Ca  ('  R.J  1*1 

00002640 

0197 

NSX*NSX.l 

00002643 

0198 

STX(NSXI=  (LR*100*JI*1000000 

00002646 

89 
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0199 

FIX( J 1-1 

00002650 

c 

LQl  KEEPS  AN  ACCOUNT  OF  COLUMNS  FOR  WHICH  FIX(J)-1 

00002655 

0200 

LQ1-LQ141 

00002660 

C 

FIXl(J)  SPECIFIES  INDEX  I  FOR  WHICH  FIX(J)-1 

00002662 

0201 

FIXK  JI-LR 

00002665 

0202 

GO  TO  1500 

00002670 

0203 

1250 

CONTINUE 

0000.’680 

0204 

1300 

CONTINUE 

00002690 

0205 

1500 

CONTINUE 

00002700 

0206 

1800 

IF  (1STEP.NE.2)  60  TO  2000 

00002710 

0207 

PRINT  1900<  K»  LQ2t  LQl 

00002720 

0208 

1900 

FORMAT  (•0‘t*K-**13.‘  L02=*,I3,  •  L01=',  13) 

00002730 

0209 

00  1930  1-ltM 

00002740 

0210 

PRINT  290t  It<CX<I,J)tJ-l.N) 

00002750 

0211 

1930 

CONTINUE 

00002770 

0212 

PRINT  297,  (FIX(JI,J-1.NI 

00002780 

0213 

2000 

CONTINUE 

00002800 

c 

A  CYCLE  EXAMINES  ALL  THE  FACILITIES. 

00002803 

c 

IF  IN  A  CYCLE,  THE  CAPACITY  RULE  RESULTS  IN  SETTING 

00002810 

c 

ACiOITIONAL  CX(I,J)  VALUES  EQUAL  TO  2,  THEN  REPEAT 

00002820 

c 

THE  CYCLE.  BUT  IF  FIXIJ)-!  FOR  ALl  J,  THEN  DO  NOT 

00002830 

c 

REPEAT  THE  CYCLE. 

00002835 

0214 

IF  ILOl.EO.N)  GO  TO  2400 

00002840 

0215 

IF  (LQ2.EQ.LR2I  GO  TO  2400 

00002845 

0216 

2200 

LR2-LQ2 

00002860 

0217 

60  TO  310 

00002870 

C 

00002880 

C 

UPDATE  VECTOR  OF  FACILITIES  FLB(K)  FOR  COMPUTING 

00002890 

c 

C(I,J1  MATRIX  C  LOWER  BOUND.  IT  HAS  VALUE  1  IF  A 

00002900 

c 

FACILITY  IS  USED,  OTHERWISE  IT  HAS  0  VALUE. 

00002910 

C21d 

2400 

DO  3000  J-1,N 

00002950 

0219 

IF  (FIX(J).EQ.O)  GO  TO  3000 

00002960 

0220 

INDI-FIXK  Jl 

00002970 

0221 

DO  2550  K-l,P 

00002990 

0222 

IF  (E(INOI,KI.EQ.O)  GO  TO  2550 

00003000 

0223 

IF  (FL8(K).E0.1)  GO  TO  2550 

00003010 

02  24 

FLB(K)=1 

00003020 

0225 

2550 

CONTINUE 

00003030 

0226 

3000 

CONTINUE 

00003060 

0227 

IF  (ISTEP.NE.2)  GO  TO  3150 

00003070 

0228 

PRINT  3100,  (FLB(K) ,K=1,P) 

00003080 

0229 

3100 

FORMAT! 'O’ ,•  (FLB  (K)  ,K»1, PI  ♦,  201 4/’l6X, 2014 1 

00003090 

C 

COMPOTE  COST  MATRIX  C(1,J)  FOR  THE  RELAXED  PROBLEM 

00003100 

0230 

3150 

DO  3400  J-1,N 

00003110 

0231 

DO  3300  I-1,M 

00003120 

02  32 

BSUM<cO.O 

00003130 

0233 

DO  3200  K«1,P 

00003140 

0234 

IF  (FLB(KI .EQ.l )  GO  TO  3200 

00003150 

02  35 

IF  (Ed  .KI.EO.O)  GO  TO  3200 

00003160 

0236 

BSUMs8SUM4|a(KI  •  (FLaAT(OII,J,K))/  FLOATISIKdll 

00003170 

0237 

3200 

CONTINUE 

00003180 

0238 

3250 

C(  I, J)-A( I , J)*BSUM 

00003190 

02  39 

3300 

CONTINUE 

00003200 

0240 

3400 

CONTINUE 

00003210 

0241 

IF  (ISTEP.NE.2)  GO  TO  3445 

00003220 

0242 

DO  3430  I-1,M 

00003230 

0243 

PRINT  3420,  I,  (C( I,J),J-1,N) 

00003250 

0244 

3420 

FORMAT  (/5X,  •C(I,J)*,5X,  •!»•,  13, 2X,  5F15.4, 

00003260 

il 
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1  6(/23X,  5F15.4M 

00003265 

02aS 

3430 

CONTINUE 

00003270 

C 

FIND  SUM  OF  MINIMUM  C(I,J>  VALUES  OVER 

EACH  J, 

00003290 

C 

I.E.t  MINSCsSUM  of  HJNC(J). 

00003300 

c 

IF  FIX(J»*lf  THEN  MINCIJ)>=CIT,J)  WHERE 

CXd,J)*l 

00003310 

0246 

3445 

MINSC=0.0 

00003320 

0247 

DO  3900  Jsl.N 

00003340 

0248 

IF  (FIX(JI.EO.O)  GO  TO  3500 

00003350 

0249 

INOI^FIXKJI 

00003360 

0250 

MINC(J)=C( INOI.OI 

00003370 

0251 

S0LX(8I>:1N0I 

00003380 

0252 

GO  TO  3850 

00003410 

0253 

3500 

LK=0 

00003430 

0254 

1  =  1 

00003440 

C 

SKIP  CdtJt  ELEMENT  IF  CX(v,J)s2  L  MOVE 

TO  NEXT  1 

00003470 

0255 

3550 

IF  (CX( l,J).Ea.2l  GO  TO  370: 

00003480 

0256 

IF  l(I-LK).EO.l)  GO  TO  3600 

00003485 

0257 

IF  (CdtJI.GE.MINCIJl)  GO  TO  3750 

00003490 

0258 

3600 

MlNC(J>=Cd.J) 

00003500 

0259 

IMlNxl 

00003510 

0260 

GO  TO  3750 

00003520 

0261 

3700 

LK=LK*1 

00003530 

0262 

3750 

1=1*1 

00003590 

0263 

3800 

IF  (I.LE.Ml  GO  TO  3550 

00003600 

02o4 

SQLX(Ji=IMIN 

00003610 

0265 

3850 

MINSC=M1NSC*MINC«J) 

00003620 

0266 

3900 

CONTINUE 

00003630 

02o7 

IF  dSTEP.NE.2l  GO  TO  3940 

00003640 

02o8 

00  2920  J=1,N 

00003650 

0269 

PRINT  3910,  J,M1NC(JJ,S0LX(J) 

00003660 

0270 

3910 

format  CO'.'J.MlNCIJl.SOLXlJl'f  15,F15.4,T6) 

00003670 

0271 

3920 

CONTINUE 

00003680 

C 

COMPUTE  FIXED  COST  FC  FOR  LC WER  BOUND 

00003710 

02  72 

3940 

FC=0 

00003720 

0273 

DO  4000  K=l,P 

00003730 

02  74 

IF  (FLB(KI.EO.O»  GO  TO  4000 

00003740 

02  75 

3950 

FC=FC*B(K) 

00003750 

0276 

/  000 

CONTINUE 

00003760 

C 

00003770 

0277 

4050 

LOWBsMINSC*FC 

00003780 

02  78 

IF  dSTEP.EQ.O)  GO  TO  4150 

00003790 

0279 

PRINT  4120,  MlNSC,  FC,  LOWS 

00003800 

0280 

4120 

FORMAT  (•O',*  MINSC,  FC,  L0W8  •,  F15.4,  115,  F15.4I 

00003810 

C 

COMPARE  LOWER  BOUND  WITH  BEST  UPPER  BOUND  STAR 

00003820 

C 

BUBS  WHICH  EQUALS  BUB/d*EPS).  IF  LOWB 

IS 

00003830 

C 

GREATER  THAN  OR  EQUAL  TO  BUBS,  THEN  BACKTRACK 

00003840 

0231 

4150 

IF  (LOWB.GE.BUBSIGO  TO  6200 

00003850 

C 

CHECK  IF  CURRENT  SOLUTION  SATISFIES  CAPACITY 

00003880 

C 

CONSTRAINTS 

00003890 

0282 

4200 

IF  dUNCAP.EC.il  GO  TO  4420 

00003900 

0283 

4210 

DO  4400  K=1,P 

00003910 

0284 

NSUMDsO 

00003920 

0285 

DO  4300  J=1,N 

00003930 

C2U6 

1X=SCLX ( J 1 

00003950 

0287 

NSUMD=NSUMO*DdX,a,  K) 

00003960 

0268 

4300 

CONTINUE 

00003970 

0289 

IF  dSTEP.NE. 21  GO  TO  4320 

00003980 

0290 

PRINT  4310,  K,NSUNO 

00003990 
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FORMAT  (‘O*.  •K,NSUMO*t21IO) 

IF(N$UMO.LE.S(K) )  GO  TO  **00 
GO  TO  5100 
CONTINUE 

«««*««*«««*«««CCMPUTE  UPPER  BOUND  UPS  IF  CAPACITY  CONSTRAINTS 
ARE  SATISFIED. 

UPB^SUM  OF  Ad.JI^FlXED  COST  FCU3  BASED  ON 
SOLUTION  VECTOR  SOLX(J) 

VECTOR  OF  FACILITIES  FOR  UPPER  BOUND  FUBIKt  HAS 
VALUES  I  OR  0  BASED  ON  FACITILV  USED  OR  OTHERWISE 

DO  4*50  Rsl,P 
FuecK i«a 
CONTINUE 
NSUMA=0 
FCUB^O 

DO  *650  J>ltN 
1X=S0LX (J) 

NSUHA>NSUM A«A ( IX  t  J i 
00  *600  K>1,P 

IFCEIIX.KI.EO.OI  GO  TO  *600 
tF(FUB<KI.EQ.ll  GO  TO  *600 
FUB(KI>1 
FCUBxFCUB«S(KI 
CONTINUE 
CONTINUE 

IF  (ISTEP.NE.2J  GO  TO  *700 
PRINT  *660f  (FUB(K) ,Kxl,P) 

FORMAT! 'O*  ,MFUB(K)  ,Kxl,P)  *,  201*/I6X,  201*) 

UPBxNSUMA*FCU8 
IF  (ISTEP.EO.OI  GO  TO  *750 
PRINT  *710t  NSUMA,  FCUBt  UPBt  BUB,  BUBS 
FORMAT! 'O' t'NSUMA,  FCUB,  UP6,  BUB,  BUBS  *,2110,  .'FIS.*) 
compare  upper  sound  WITH  BEST  UPPER  BOUND 
IF  UPB  IS  LESS  THAN  BUB,  SET  IT  AS  BUB  AND 
NOTE  THE  SOLUTION 
IF  lUPB.GE.BUB)  GO  TO  5100 
BUBxUPB 

BUBS«  BUB/  !1.0*EPS) 

IBNOOxNOO 

PRINT  *780,  IBNOO,  BUB,  BUBS 

FORMAT  !'0',  MBNOO,  BUB,  BUBS',  110,  2F15.*) 

□0  *800  Jxl,N 
BSOLX! JIxSOLX! J) 

DO  *850  K>1,P 
8S0LY!K|xFUB!K) 

•••««*««***«««COMPARE  LOWB  WITH  BUBS.  IF  LOWB  IS  GREATER 
THAN  OR  EQUAL  TO  BUBS,  THEN  BACKTRACK 
IF  !LOWB.GE.BUBS)GO  TO  6200 

•««««***«**««*IF  FIX!U)  VALUES  ARE  1  FOR  EACH  J,  THEN  BACKTRACK 
IF  !L01.E0.N)  go  to  6200 

THE  BOUNDING  RULE*****'******''************ 

IF  THE  DIFFERENCE  BETWEEN  C!I,J)  AND  H1NC!J)  IS 
GREATER  Than  THE  DIFFERENCE  BETWEEN  BUBS  AND 
LOWB,  THEN  CK!1,J)*2 

••••••••••**P**PPLV  BRANCHING  RULE  AND  FIND  SRI,  THE  NEXT 

VARlAeLE  FOR  LEFT  BRANCHING. 

FIND  NMINCIJ),  THE  NEXT  HIGHER  VALUE  THAN  MINC!J) 
ANO  OIFBRiJ),  THE  DIFFERENCE  BETWEEN  THEM. 


0000*000 
0000*010 
0000*020 
0000*030 
0000*0*0 
0000*050 
0000*060 
0000*070 
0000*080 
0000*090 
0000*100 
0000*110 
0000*120 
00004130 
0000*1*0 
OOC  *150 
0000*170 
0000*180 
0000*190 
0000*200 
00004210 
00004220 
00004230 
000042*0 
0000*250 
00004260 
0000*270 
0000*280 
0000*290 
00004300 
0000*310 
0000*320 
00004330 
000043*0 
0000*350 
0000*360 
0000*370 
00004380 
0000*385 
0000*386 
0000*388 
0000*390 
00004*00 
0000**10 
0000**20 
0000*430 
0000***0 
00004*50 
0000*480 
0000*500 
0000*510 
00004515 
0000*520 
0000*525 
0000*530 
0000*5*0 
0000*550 
0000*555 
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huKIRAN 

0329 

0330 

0331 

0332 

0333 

0334 

0335 

0336 

0337 

0338 

0339 

0340 

0341 

0342 

0343 

0344 

0345 

0346 

0347 

0343 

0349 

0350 

0351 

0352 

0353 

0354 

0355 

0356 

0357 

0358 

0359 

0360 

0361 

0362 

0363 

0364 

0365 

0366 

0367 

0368 

0369 
0370 
0371 
0372 
0373 
03  74 
0375 
0376 
0377 
0378 
03  79 
0330 
0331 
03  32 
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5300 

5320 

5330 


DBOUNO>eUBS-LOMB 
5200  00  5250  J*lfN 

NH1NC( J>s0.0 
DIEBRIJ)=0.0 
5250  CONTINUE 

00  5600  J-lfN 

C  SKIP  TO  NEXT  J  IF  FIXJJ)*! 

IF  (FIXIJI.EQ*!!  GO  TO  5600 
uK«0 

C  SKIP  Cd.JI  IF  CXUtJI«=2  C  /<QVE  TO  NEXT  I 

IF  (CX(I,J).E0.2)  GO  TO  5350 
IF  (I.EQ.SOLXIJ) )  GO  TO  5350 
IF  {JCdf  J»-M1NC(JI  l.GT.DBOUNO)  GO  TO  5330 
IF  U1-LK).£0.1)  GO  TO  5320 
IF  K(ItJ).GE.NMINC(jn  GO  TO  5400 
NMINC(JI»C(ItJ) 

GO  TO  5400 
CX(I .J)-2 
NSX=NSX*1 

STX<NSX>=-(1*1004JJ-1000000 
KT2(J»*KT2(J»*1 

lF(KT2(Jt.LT.(H-lH  GO  TO  5350 
1NDI=S0LX( J) 

CXIINDI,J>=1 
N'.X=NSX*1 

STX(NSX|s  (IN0I*100»J>*1000000 
FIX(J»=1 
L01-LQ141 
FIXK  J>>1N0I 
GO  TO  5600 
LK*LK4l 
I»l4l 

IF(I.LE.M)  GO  TO  5300 
01FBR(J)=NHINC(J)-HINCIJ) 

CONTINUE 

IF  (ISTEP.Ne.2l  GO  TO  5650 
00  5620  1»1|H 

PRINT  290,  I,ICX(I,J»,J*l,N» 

CONTINUE 

PRINT  297,  <FIX(J)  ,J»1,NI 

IF  FIX<JF=I  FOR  *LL  J,  THEN  eACKIRACK. 

IF  (LQl.EQ.NI  GO  TO  6200  _ 

FIND  HAXOIF,  THE  MAXIMUM  DIF-ERENCE  DIFBRIJ) 

LF=0 

00  5800  J=1,N 

IF  (FIXI Ji .EO.l I  GO  TO  5690 
IF  HJ-LF».EO.H  GO  TO  5660 
IF  (OIFBRC JI .LT.MAXDIFI  GO  TO  5800 
5660  MAXOIFsOIF8R( JI 


5350 

5400 

5500 

5600 


5620 


5650 


LJ-J 

GO  TO  5800 
5690  LF*LF*1 
5800  CONTINUE 

IF  !IST£P.NE.2I  GO  TO  5840 
DO  5820  J>1,N 

IF  IFIXiJi  .EO.l  I  GO  TO  5820 

PRINT  5810,  J,  NMINCCJI,  MINCIJI,  DIFBRIJI 


00004568 

00004570 

00004580 

00004590 

00004600 

00004610 

00004620 

00004630 

00004640 

00004650 

00004670 

00004680 

00004690 

00004700 

00004710' 

00004720 

00004730 

00004735 

00004740 

00004742 

00004745 

00004747 

00004750 

00004752 

00004755 

00004758 

00004760 

00004762 

00004764 

00004766 

00004768 

00004770 

00004775 

00004780 

00004785 

00004790 

00004795 

00004820 

00004830 

00004850 

00004860 

00004880 

00004890 

00004900 

00004905 

00004910 

00004915 

00004920 

00004925 

00004930 

0000493b 

00004940 

00004943 

00004946 

00004950 

00004953 

00004956 

00004960 
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0383 

5810 

FORMAT  (•0't*J.NMINC(Ji<MtNCIJIt01FBR|J)',  15,3F15.4) 

00004963 

0384 

5820 

CONTINUE 

00004966 

C 

00004970 

0383 

5840 

DO  5900  J»1,N 

00004980 

0386 

IF  (J.NE.LJ)  GO  TO  5900 

00004990 

0387 

5850 

BR1=S0LXIJI«100aJ 

00005000 

0388 

IF  IISTEP.EQ.O)  GO  TO  6020 

00005010 

0339 

PRINT  5880,  BRl 

00005020 

0390 

5880 

FORMAT! 'O', •  BRl*. 1101 

00005030 

0391 

GO  TO  6020 

00005040 

0392 

5900 

CONTINUE 

00005050 

C 

00005060 

C 

NSX  REPRESENTS  THE  NUMBER  OF  VARIABLES  IN  STXMNS) 

00005070 

0393 

6020 

NSX=NSX*1 

00005090 

0394 

6040 

STXINSX)>BR1 

00005100 

0395 

IF  (ISTEP.NE.2t  GO  TO  6100 

00005150 

0396 

PRINT  6088,  (STXdNSI,  INS  =  1,NSXI 

00005160 

0397 

o088 

FORMATCO*,*  STXdNSI*,  10110,  1221/,  lOX.lOIlOII 

00005170 

C 

00005220 

C 

CAPACITY  RULE 

00005230 

0398 

6100 

NOO=NQO«1 

00005240 

0399 

6110 

IF  (ET.EO.O.OI  GO  TO  6150 

00005242 

0400 

IFdNSET.EO.il  GO  TO  6147 

00005244 

0401 

IF  (INET.EO.ll  GO  TO  6150 

00005246 

0402 

CALL  TIMETdNTI 

00005248 

0403 

ELTN=dNT-lI0l*26.04E-6 

00005250 

0404 

IF  IELTN.lt. ETI  GO  TO  6150 

00005253 

0405 

6120 

PRINT  6125,  NOD,  ELTN,  BUB,  BUBS,  IBNOO 

00005256 

0406 

6125 

FORMAT  (*0*,  *WAS  AT  NODE*, 16,  *  AT  ELAPSED  TIME  =*,  F10.4, 

00005260 

1 

*  SECONDS. ' ./IX,  •  8U6**,F15.4,  *  BUBS** ,F15.4, 

00005263 

2 

*  AT  NODE**, 171 

00005266 

0407 

IBUB*8UB 

00005267 

0408 

IF  (IBUB. EQ. 9999999 1  GO  TO  6l46 

00005268 

0409 

6130 

PRINT  6135,  (BSOLXIJI . J=1,N» 

00005270 

0410 

6135 

F0RMAT«*0',  *SOLUTION  CORRESPONDING  TO  BUB  IS*,  //IX, 

00005273 

1 

*(BSOLX(JI,  J*1,NI*, 1018, 3I/18X, 10181) 

00005276 

0411 

6140 

PRINT  6145,  «BS0LY(K1,  K=1,P) 

00005280 

0412 

6145 

FQRMAT|/lX,*(aSOLY(KI,  K =1 ,P 1  * , 10 18 ,2 ( /I 8X ,  lOlBI) 

00005290 

0413 

6146 

INET*! 

00005292 

0414 

INIS=1STEP 

00005294 

0415 

INSET=1 

00005296 

0416 

1STEP=2 

00005298 

0417 

GO  TO  6150 

00005300 

0418 

6147 

ISTEP=INIS 

00005302 

0419 

INSET=0 

00005304 

0420 

6150 

GO  TO  272 

00005306 

C 

00005308 

0421 

6200 

IF  (NSX.EQ.O)  GO  TO  8100 

00005310 

0422 

6220 

IF  (  lABSISTXlNSXl 1 .GT. 1000000)  GO  TO  6500 

00005320 

0423 

62  50 

BRO=STX(NSX) 

00005330 

0424 

6270 

STXINSX )*-BR0-l0C0000 

00005340 

0425 

IF  dSTEP.tO.OI  GO  TO  6308 

00005390 

0426 

PRINT  6305,  BRO 

00005400 

0427 

6305 

FQRMAT(*0*,  'BRO  *,110) 

00005410 

0428 

6308 

IF  dSTEP.NE.2l  GO  TO  6330 

OOOOSwiO 

0429 

PRINT  6088,  (STXdNSI,  INS*I,NSX) 

00005430 

C 

00005490 

C 

CAPACITY  RULE 

00005500 

94 
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0430 

6330 

NCD=NaO*l 

00005510 

0431 

6410 

If  (er.EO.o.o)  GO  ro  645o 

00005512 

1 

0432 

IF  (INSET. EQ.l)  GO  TO  6445 

00005516 

0433 

IF  (INET.EC.l)  GO  TO  6450 

00005518 

0434 

CALL  TIHETdNTI 

00005520 

( 

04  35 

ELTN«llNT-IT0)*26.04£-6 

00005523 

0436 

IF  (ELTN.LT.tT)  GO  TO  6450 

00005526 

0437 

6420 

PRINT  6125«  NOO.  ELTNt  BUB,  BUBS,  IBNOO 

00005528 

0 

0438 

IBUBxBUB 

00005530 

0439 

IF  (1BUB.E0.9V99999)  GO  TO  6442 

00005532 

0440 

6430 

PRINT  6135,  (BSOLXt J),J=1,N) 

00005533 

& 

0441 

6440 

PRINT  6145,  (BSOLY(KI,  K>sl,P) 

00005536 

0442 

6442 

INET=1 

00005538 

0443 

INIS^ISTEP 

00005540 

o 

0444 

INSET=1 

00005542 

0445 

JSTEP*2 

00005544 

0446 

GO  TO  6450 

00005546 

o 

0447 

6445 

ISTEP=INIS 

00005548 

0448 

lN5eT=0 

00005550 

0449 

6450 

GO  TO  210 

00005552 

c 

0450 

6500 

IF  (  STX(NSX) .GT. 10000001  GO  TO  6520 

00005555 

0451 

LX=-STX<NSX (-1000000 

00005560 

0452 

IX=LX/100 

00005570 

0453 

JX»LX-IX*100 

00005580 

0454 

CX(IX,JX)=0 

00005590 

0455 

Kr2(JXI»KT2(JX(-l 

00005595 

r 

0456 

CO  TO  6550 

00005600 

0457 

6520 

LX*  STX<NSX)-1000000 

00005610 

0458 

IX=LX/100 

00005620 

• .  . 

0459 

JX=LX-IX*100 

00005630 

0460 

CX(1X,JX)*0 

00005640 

0461 

F.X( JX)=0 

00005650 

■J 

0462 

LQ1*L01-1 

00005660 

0463 

6550 

NSX=NSX-1 

00005690 

0464 

GO  TO  6200 

00005700 

C 

00005730 

0465 

8100 

16ue*6UB 

00005740 

0466 

CALL  TIMETdTlI 

00005750 

O 

0467 

ELTl*dTl-IT0)*26.04E-6 

00005760 

0468 

PRINT  8105,  ELTl 

00005770 

0469 

8105 

FORMAT  CO*, ///IX,  'ELAPSED  TIME  IN  SECONDS-',  F15. 

8) 

00005780 

0 

0470 

PRINT  8120,  NOO 

00005790 

0471 

8120 

format  (  'O', 'TOTAL  NUMBER  OF  NODES  EXPLORED  »*,I3I 

00005800 

0472 

IF  dBUB.EQ. 9999999)  GO  TO  8350 

00005810 

c 

0473 

8130 

PRINT  8150 

00005820 

0474 

8150 

format  ('O',  'NOTE:  1.  FOLLOWING  X(I,J) 

variables 

SHOW  DESIGN', 

00005830 

1  •  I  TO  WHICH  Acrivirv  j  is  assigned  for  J»l 

TO  N. ', 

00005840 

9 

2  /7X,  '2.  IF  EPSILON  EPS  WAS  ASSIGNED  A  POSITIVE', 

00005850 

3  •  (NON-ZERO)  value,  THE  SOLUTION 

MAY  BE  SUBOPTIMAL.',/) 

00005860 

0475 

8180 

PRINT  8200,  (8S0LX(Jt,J*l,N) 

00005870 

t 

04  76 

8200 

FORMAT! *0' ,155,  'OPTIMAL  SOLUT ION' ,/ 1 X , T 

J5, 

00005B80 

2  3(/,21X,10iai) 

00005900 

0477 

8210 

PRINT  8250,  (8S0LY(M,  K*I,P) 

000059)0 

04  78 

6250 

FORMAT  ('O',  'Y(K)  VALUES:',  8X ,  1018 

,  2(/,21X, 

loiai ) 

00005920 

0479 

6280 

PRINT  8300, IBU8 

00005930 

( 

04d0 

8300 

format  ('O',  'OPTIMAL  VALUE  OF  OBJECTIVE 

FUNCTION: 

',  115///) 

00005940 

0481 

GO  TO  8500 

00005950 

I 
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0Hd2  8350  PRINT  8400 

0^83  8400  FORMAT  I'O*#  •  PROBLEM  DOES  NOT  HAVE  A  FEASIBLE  SOLUTION*, 

1  /IX,  •  BECAUSE  THE  CAPACITY  CONiTRAINTS  CANNOT* 

2  /IX,  •  BE  SATISFIED. *,/! 

0484  8500  PRINT  8550 

0485  8550  FORMAT  ( *0  * , **P***NORMAL  END  OF  JOB***PP*,/) 

U486  8600  STOP 

0487  END 


00005960 

00005970 

00005980 

00005990 

00006000 

00006010 

00006020 

00006030 


APPENDIX  B 


DETAILED  PRINTOUT  FOR  A  TEST  PROBLEM 
(TEST  PROBLEM  WITH  m=5,  n=4,  AND  p=&> 
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